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Abstract 

We analyze isomonodromic deformations of rational connections on the Ricmann sphere with Fuchsian 
and irregular singularities. The Fuchsian singularities arc allowed to be of arbitrary resonant index; the 
irregular singularities arc also allowed to be resonant in the sense that the leading coefficient matrix at 
each singularity may have arbitrary Jordan canonical form, with a genericity condition on the Lidskii 
submatrix of the subleading term. We also give the relevant notion of isomonodromic tau function 
extending the one of non-resonant deformations introduced by Miwa-Jimbo-Ueno. The tau function is 
expressed purely in terms of spectral invariants of the matrix of the connection. 
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1 Introduction 



The history of isomonodromic deformations dates back to Schlcsingcr [17] who studied rational connec- 
tions of size r x r with simple poles (Fuchsian singularities) 



and defined a system of differential equations describing the dependence of the residue matrices Ai on the 
position of the poles 7, under the condition that the monodromy representation induced by the kernel of 
the connection is independent of the position of the poles. 

Later on the authors of [11] described an extended system of equations where irregular singularities 
were considered; 



in this case the differential equations are with respect to certain exponents in the formal asymptotic 
behavior of the kernel solution and the "monodromic" data are extended to include the parameters 
appearing in Stokes' phenomenon. 

In both settings the study was limited to non-resonant connections. The notion of resonance for a 
rational connection is as follows; for a Fuchsian singularity we say that the connection is resonant if the 
residual spectrum at a singularity (the spectrum of the residue matrix) contains eigenvalues differing by 
a nonzero integer. An irregular singularity is called resonant if the eigenvalues of the leading coefficient 
are repeated. 

In the past years significant classes of examples have appeared in which one is forced to consider 
resonant isomonodromic deformations. Probably one of the most relevant is in the analysis of quantum 
cohomology using the notion of Frobenius manifolds [8] , where a Fuchsian resonant singularity appears 
in the relevant isomonodromic deformation. 

Lately, the analysis of the Ricmann-Hilbert problem associated to biorthogonal polynomials for the 
multi-matrix models has lead to a system with an irregular resonant singularity [3, 4]. 

The general structure of an isomonodromic deformation of a resonant Fuchsian connection was ad- 
dressed in [7] whereas it seems that no attempt is being made in the literature to address isomonodromic 
deformations of irregular resonant singularities in general. It is the purpose of this paper to analyze these 
issues. The completely general classification of Fuchsian resonant singularities and their isomonodromic 
deformations are essentially contained in [7] and we are rephrasing it in Section 4 for the reader's conve- 
nience. The main feature which distinguishes these deformations from the usual Schlesinger equations is 
that in a monodromy-preserving deformation of a resonant Fuchsian connection the deformation matrices 
may in fact have higher order poles, of degree at most equal to the maximal integral difference between 
two residual eigenvalues at each pole (the resonance index) ; these higher order poles -however- are the 
result of the bigger "local" gauge freedom that arises due to the resonant character of the singularity and 
can in fact always be gauged to zero by a rational gauge equivalence (without changing the position of 
the poles and the Fuchsian character of the connection, see Sect. 4). 

The situation is not dissimilar from the nonresonant case, in which -however- the gauge freedom is 
restricted to a point where the only arbitrariness is global constant gauge transformations. 

The case of isomonodromic deformations of irregular resonant singularities (in the generalized sense 
of [11]) is quite unexplored, mainly because of the difficulty in analyzing the normal asymptotic form 
near any such singularity. 

The class of resonant singularities which we analyze here may well be considered "minimally" resonant 
in the sense to follow: we will consider rational connections A{x) such that near an irregular singularity 
x = 7 the leading coefficient matrix A r _ 7 may have an arbitrary Jordan canonical form. However we 
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impose a (completely explicit) genericity assumption on the second-leading coefficient matrix as explained 
in Section 2. 

The meaning of our genericity assumption has a clear interpretation in terms of the spectral curve 
of the connection A(x), i.e. the algebraic curve satisfied by the eigenvalues y(x) of A{x). It ensures 
that the branching structure of the (desingularized) spectral curve above x = 7 is in agreement with the 
dimension of the Jordan cells of the leading coefficient matrix. 

This allows us to obtain a canonical form for the asymptotic behavior of a kernel solution ty' = 
(Thm. 3.1), which is the crucial tool in order to construct isomonodromic deformations and prove their 
compatibility. 

Tau Function. Another main point of our construction is that we define -following [11]- the notion of 
"isomonodromic tau function" for the deformations we have defined. Here our approach differs radically 
from that in [11] inasmuch as our definition of tau function is not obtained in terms of the formal 
asymptotic data. We rather use the spectral curve itself, thus showing explicitly the spectral nature of 
the tau function. The approach is along the same lines of [6], where it was shown that Miwa-Jimbo-Ueno's 
tau function is in fact a spectral invariant expression. 

Let us briefly comment on the necessity and interest in a definition of tau function in this context; it 
has been shown repeatedly [2, 5, 9] that JMU's tau function coincides with Toplitz/Hankcl determinants 
of moments of measures. This allows to establish a direct connection between the partition function of 
certain matrix models and the tau function of a (naturally) associated isomonodromic deformation. 

In the context of multimatrix models we have many of the same features, namely a partition function 
and isomonodromic deformations of a certain ODE [3, 4] and it is natural to imagine an analogous 
relationship with an isomonodromic tau-function. The main obstacle is the absence of a general definition 
of tau function for resonant irregular singular ODEs. 

Although the resonance of the ODE appearing in the two-matrix model is not quite of the class 
which we analyze in this paper, on the other hand it has many of the same features. In particular the 
tools developed here are sufficient to analyze that situation (although in a ad hoc way) . This analysis is 
contained in a separate paper. 

2 Preamble: perturbations of spectra and Lidskii coefficients 

In this section we introduce some necessary notations and definitions that will be used in the analysis of 
the formal asymptotics of singular ODEs. The philosophy inspiring these considerations is that solving 
a singular ODE by formal series is -to high order- the "same" as finding perturbative eigenvectors of an 
analytic (formal) perturbation. We recall some relevant facts about analytic (or formal) perturbations of 
the spectrum of matrices. From our point of view the issue of convergence of the series that will appear 
is irrelevant since that in our applications only a finite number of the coefficients will appear; for this 
reason we will limit the discussion to the formal aspects of the problem, with the understanding that 
under mild additional assumptions the considerations to come could be set in an analytic framework. 
By perturbation we mean a (formal) power series in a small parameter e of the form 



One of the main question in perturbation theory is to understand the behavior of the spectrum of M(e) 
and its relation to the "unperturbed" spectrum of M(0). A related question is that of describing the 
perturbation of the corresponding eigenvectors. 

The generic case is very well understood and studied and corresponds to the case where all unperturbed 
eigenvalues are distinct and simple (i.e. with algebraic multiplicity one). In this case it is not hard to 
show that each perturbed eigenvalue admits a power series expansion (possibly formal) [12] 



M(e) := M + eMi + C(e 2 ) . 



(2-1) 



A J (e)=A,(0)+O(e). 



(2-2) 
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Complications arise when the unperturbed spectrum consists of eigenvalues with algebraic multiplicity 
higher than one, namely when the Jordan canonical form of Mo is allowed to be the most general. 
Without loss of generality we may assume that Mo is in its Jordan canonical form. Adopting Arnold's [1] 
notation we denote a matrix in Jordan canonical form by the product of the determinants of its blocks. 
For example a 3 a 2 {3 denotes a matrix with two Jordan blocks with eigenvalue a of dimension 3x3 and 
2x2, and another Jordan block of size 4x4 with eigenvalue (3. In this example then the algebraic 
multiplicity of a is 3 + 2 = 5 but the geometrical multiplicity (the rank of the eigenspace) is 2. 

As a general rule a multiple eigenvalue like the a in the example, will split under perturbation in 5 
distinct eigenvalues; in fact under certain gencricity assumption on the first jet of the perturbation (i.e. 
Mi) the splitting that will occur is in a triplet and a doublet according to the size of the two Jordan 
blocks. The eigenvalues of the triplet will be Puiseux series in ea and will be cyclically permuted after 
a loop around the origin in the e plane; similarly the eigenvalues of the doublet will be Puiseux series in 
e 1 / 2 enjoying a similar cyclicity. 

In [14] were studied sufficient conditions for the splitting of an eigenvalue to occur in cyclic fc-tuplcts 
according to the sizes of the Jordan blocks of the unperturbed matrix Mo- More recently the approach 
of Lidskii has been refined (see [10] and references therein) to handle the cases in which the splitting of 
eigenvalues does not necessarily respect the Jordan decomposition of Mo . For the purposes of the present 
paper we restrict ourselves to the "generic" stratum within this class. 

In order to describe this theory we start by observing that we can limit ourselves to the case where 
Mo has a single eigenvalue which we can set to zero by shifting with the identity matrix. Indeed blocks 
with distinct eigenvalues will have eigenvalues which will not "mix" under a small perturbation. 
Suppose thus that Mo = ni n2 ■ ■ ■ ,lK ; we arrange these blocks in weakly decreasing order 

ni > n>2 > ■ ■ ■ > tik- (2-3) 

We next partition the first jet of the perturbation (Mi) in K x K blocks according to the same block 
decomposition of Mo; from each of these blocks of Mi we extract the lower-left entry and form a K x K 
matrix; we will call such matrix L(Mi)/ ni) „ 2j ... iraK } the Lidskii matrix of Mi subordinated to the block 
decomposition of Mq. 
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Mi = 



L21 



L 



L 



22 



13 



-23 



Li 



L 



24 



(2-5) 



^31 
-41 



-32 



^33 
-43 



L34 
L44 



Let £\ be the number of blocks of the same size as n\ : the next block of strictly smaller dimension will 
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be then the £\ + 1. Let £ 2 be the number of blocks of the same size as the n^+i, and so on and so forth. 
At the end of this procedure the K diagonal blocks are grouped together according to the dimensions 

n-t,n 2 , ...,ni 1 ,m 1+ i,...,ni 1+ i 3 ,... . (2-6) 

V v ' V V ' 

This grouping induces a partitioning of the Lidskii matrix L into blocks; for the example in the figure 
above we have (£1 = 1,£ 2 = 2, £3 = 1) 





' in 


L12 


L13 


Ll4 


L = 


L21 
L31 


L22 
L32 


L23 
L33 


L24 
L34 




L41 


L42 


L43 


L44 



(2-7) 



We now consider the principal block-submatrix of L according to this decomposition: namely the first 
principal block submatrix is of size £±x£i, the next is of size (£i+£2) 2 , etc. For each of these submatrices we 
construct the pseudo-characteristic polynomial, namely the determinant of the submatrix minus A times 
the projector onto the lower right sub-matrix. At each step we have a submatrix of size (£1 + . . . + £j) 2 
and the corresponding pseudo charpoly is of degree £j . 

Definition 2.1 The roots of these polynomials will be called the Lidskii pseudovalues. 

We will say that they are generic if none of them is zero and the discriminant of each pseudocharpoly 
is nonzero. 

We think that the description of the procedure is sufficiently involved to require an example: contin- 
uing with the above one, the Lidskii pseudo-charpoly's are 



Pi 


= det[L 


11 - A] 






(2-8) 








L n 


L\2 


L\3 






Pi 


= det 




L21 
L31 


L22 — A L23 
L32 L33 — A 




(2-9) 








Ln 


L12 L13 


L\4 






p 3 


= det 


L21 
L31 


L22 L23 
L32 L33 


L 2 4 

L34 


(2-10) 






L41 


L42 L43 


L44 — A 





Let us denote Aj iP , p = 1, ...,£j the Lidskii pseudovalues: then the eigenvalue of Mq splits into 

1 

K cyclic multiplets, the j-th of which is expandable as a Puiseux series in e™j and with the leading 
coefficients of these series given by the Lidskii pseudovalues. 
In the above example: 

we have a rii = 4-tuplet of cyclic eigenvalues Ai(e) = Ai^e 1 / 4 + ....; 

we have £2 = 2 ni = 713 = 2-tuplets (doublets) with expansion A2,i(e) = \i,\\ft + ■ ■ ■ and A2,2(e) = 
\2,2\fe + ■ ■ ■; 

we have £3 = 1 eigenvalue with the form A34 = A 3: ie + C(e 2 ) . 

It should be clear now that the genericity condition is such to ensure that the (germ of) the spectral 
curve at e = can be minimally resolved. 

For later purposes we now investigate a bit closer the structure of the perturbative expansion in the 
case of a single Jordan block with nonvanishing Lidskii coefficient. 

Proposition 2.1 Let M(e) = JV + YlJLi ^Mj be a (formal) perturbation of the single nilpotent Jordan 
block M of size n x n. Suppose that the Lidskii pscudovaluc is nonzero, namely Ai" := (Mi)„i 7^ 0. 
Then: 
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1. There exists a similarity transformation constant in e which transforms the problem in the following 
perturbation problem 



M(e) = 



Ai A 2 
Ai 




A n -2 A„_i 

A n -2 



A 2 
Ai 




7k 

Ai * 



+ 0{e 2 ) (2-11) 



... 

Such similarity is unique up to the centralizer of M . 

2. The coefficients X 3 ■, j = 1, . . . n — 1 are the first coefficients in the Puiseux expansion of one of 
the n cyclically permuted eigenvalues in powers of £ = e». Furthermore they depend only on the 
coefficients of Mi, rationally in Ai and polynomially in the other coefficients. 

Proof. 

Let us define 



G:=diag(0,l,2,...,n-1) 

Aj := (Mi) nl ». 
We first "shear" the perturbation problem 

s = (\ 1 o- g m(XiO g = 



(2-12) 
(2-13) 
(2-14) 



e 



where only the marked entries in the expansion above can be nonzero and the matrix C is the cyclic 
permutation 

1 

1 



C 



1 



(2-15) 



Note that the sheared perturbation problem has the following periodicity and structure 

oo 

s(o = EW 

3=1 

Si := Ail, Sj = diagonal matrices , j > 2 



(2-16) 

(2-17) 
(2-18) 
(2-19) 
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Moreover the first S2, • ■ ■ , s„ are the stars of the expansion above and involve only the corresponding 
entries of Mi. In particular the first n — j coefficients of Sj are zeroes (for j < n). 

The reason for the shearing is that the leading coefficient of this new perturbation problem has 
nondegenerate spectrum \xU) J , j = 0, . . . n and hence can be dealt with usual techniques. In fact, let us 
introduce the following eigenvector matrix W , of C 

cw = wo.- 1 , Wij = w-c*- 1 )^- 1 ) . (2-20) 

We claim that we can find a unique perturbative eigenvalue matrix of the following form 

oo 

P(0- Y^^C'pAV (2-21) 

j=o 

Po = 1 ,Pj = diagonal traceless, (2-22) 
which diagonalizes S(£) with cyclically permuting eigenvalues 

Vj (0= y&Q (2-23) 

OO 

I/(0:=E^j (2-24) 

%+1 (0=%K)- (2-25) 
In order to show this we have to solve the following formal power series identity 

s(om = p(0H0 , HO ■■= y^-'O , (2-26) 

where the unknowns are the diagonal traceless matrices pj,j > 2 and the scalars Xj , j > 2. Plugging 
the Ansatz and comparing the coefficient of the power £ K +i we obtain the recurrence relation 

K K 

sfi K p K - 3 W = ^ A„(' K •*' j) ; K ,U<> •' . (2-27) 
Solving for pk-i and using Wfl^i = C^W we obtain 

K 



Ai (pk-i - c-V-ic) = J2 (ifi-hx-iC* - C- K s 3 C K p K -^ 



(2-28) 



This recurrence relation admits a unique solution: indeed both sides are diagonal matrices and we must 
guarantee that the RHS is traceless (since the LHS is, whether pk-\ is traceless or not). This condition 
fixes Xk uniquely to be (recalling that po = 1 and all the other p's are traceless) 



K-2 , , 
X K = Tr(s K ) + ]T Tt\C~ K 8jC K p K - S J 

3=2 ^ ' 



K-2 

(2-29) 



Moreover this shows that Ai, . . . , A„ depend only on s±, . . . s n and hence only on the entries of Mi in a 
polynomial way w.r.t all entries except (Mi)„i, w.r.t which they depend rationally. 
We now revert to the original problem by inverting the shearing 

M(e) = (X^fSiOiX^ ■ (2-30) 
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The eigenvector/eigenvalue problem for S(£) is conjugated as well into the following equation which 
involves only integer powers of e 

M(e)P(e) = P{e)D(e) 

where 



P(e) 



OO 1 



-n Al J 



(2-31) 
(2-32) 



£>(e) := (A 1 C) G ^A(e)^ 1 (A 1 0- G =E-^-^"(6) 

7 — 1 ^ ' 



(2-33) 



Here we have set 



^(e) := A 1 e(AiO G C(A 1 C) 



-G 



(Ai)"e 



(Mi)„ie 



(2-34) 



Note now that the matrix P is nonsingular and (formally) analytic in e; the leading coefficient is an 
invertible constant upper triangular matrix of the form 



.7 Pj 

(Ai) 



3=1 



+ C(e) =T + C?(e). 



(2-35) 



The similarity we are looking for is (Ai) T(Ai) , where T is defined here above. In fact we have the 
leading coefficient 

(Ai) G (T- 1 M(e)T) (Ax)- G - (Ai) G (P^MPjo (Ai)- G = (Ai) G D (Ai)" G 

An-l 



(2-36) 



(Ai) C 



1 
1 



A? Af 
A • 



1 



3 
i 





(AJ- G 



Ai A 2 

Ai 



... 



An-2 A„_i 
A n - 2 







A 2 
Ai 




(2-37) 



As for the subleading coefficient it is easy to see that the (n, 1) entry is transformed to Ai, thus completing 
the proof. Q.E.D. 

Before proceeding we observe that once the perturbation matrix is in the form advocated in Prop. 
2.1 the "pseudo-eigenvalue" problem can be recast as in the following corollary 

Corollary 2.1 Suppose we have a formal perturbation problem in the form guaranteed by Prop. 2.1; 
then there exists a pseudo- eigenvector matrix of the form 



P(e) = l + 0(e) 



such that 



P- 1 MP = J2^j 

3=1 



(2-38) 



(2-39) 



8 



3 Formal Asymptotics 

The cornerstone of our analysis of resonant irregular singularities in our class is the following theorem, 
which displays the normal formal asymptotic form of a kernel solution = A*$> in a sectorial neighbor- 
hood of the singularity. 

Theorem 3.1 (Main Theorem) Let A(x) be a M x M matrix with coefficients polynomial in x of 
degree at most r — 1 (r > 1), and formal series in x~ l 



j<r 



(3-1) 



We assume that the leading coefficient is in Jordan canonical form with the following elementary block 
structure 



(Ai ni )".(A. n ') 



(3-2) 



where the eigenvalues are not assumed to be distinct. Under the genericity assumption of Def. 2.1 on the 
Lidskii pseudovalues there exists a formal gauge Y analytically invertible at oo of the form 



Y(x) =Y + Y, > det(Fo) + , 

such that the gauge transformed connection 

D{x) := Y~ X AY - Y~ l Y' = diag(Z?i, ...,D S ) 



(3-3) 



(3-4) 



is in block diagonal form according to the minimal block decomposition of A r and where a block of size 
nj corresponding to an eigenvalue \j has the form 



rrij-L 

3 j=o 



9i 

X 



HjXj 



J=0 



Gj := diag(0, 1, . . 


.,n 3 


-1) 




' 





x 




1 








n 3 (x) = 


1 


























1 



G Mat (n j, nj, C) 



The formal gauge Y is uniquely determined. 

The formal solution of the system can then be written in the form 



*form = Y(x)-V bare (x) 



bare 



diag cxp E 



H{ J x- 1 ^ , . . . , exp '-TV '] x 



(3-5) 

(3-6) 
(3-7) 

(3-8) 



(3-9) 
(3-10) 



The proof of this theorem is accomplished in various steps contained in Propositions 3.2 and 3.3, which 
consist of a refinement of the standard splitting lemma, followed by Prop. 3.4 together with Prop. 3.5, 
which specify the canonical form of the connection and its solution. 
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Remark 3.1 We are here dealing only with singularities at x = oo but an analogous statement can be made 
for (formal) series in any local parameter (x — 7) by a linear fractional transformation of the variable x in Thm. 
3.1. 

In order to prove Thm. 3.1 we show that by a formal gauge transformation we can split the connection 
A into block diagonal form according to the block diagonal structure of the Jordan canonical form of its 
leading coefficient (under the assumed genericity assumption for the subleading coefficient). 

As a preliminary step we apply the standard splitting lemma [19] which guarantees the existence of 
a formal gauge Y such that the transformed connection 

BW(x) =y~ 1 ay -y- 1 y' , 

is block diagonal and each block has a leading term with only one eigenvalue 



(i) 





V 







B, 



(i) 



... \ 



0' B$J 



K is the number of groups of coinciding eigenvalues in (3-2), the size of each block B± is the rank of the 
nilspace of the matrix A r with the same eigenvalue and different block correspond to distinct eigenvalues 
of A r . 

This can be shown by using the following well-known result inductively. 

Proposition 3.1 Suppose that the eigenvalues of A r in the series of the form (3-1) has two groups of 
eigenvalues Ai, . . . , X p and A p +i, . . . , Xk such that X{ =/= Xj for i < p and j > p. Then A r is similar to a 
block diagonal matrix 



A t 






A 2 



where Ax, A2 are of dimensions p x p and (k — p) X (k — p) respectively. 
Furthermore there exists a formal power series Y such that 



B {1 \x)=Y- x AY -Y- X Y' 



has the block diagonal form 







(i) 



B 



(i) 



where b[ , B^ are square matrices of dimensions mi x mi and m% x 1112 respectively. 
The proof of this proposition can be found for example in [19] or [18]. 

In view of the above splitting lemma we can restrict ourselves to the case where the Jordan form of 
A r has only one eigenvalue, that is 



A r = (A ni )---(A" S ) 



(3-11) 



and, without lost of generality, we can assume that A = by performing a scalar gauge transformation 
e~~ x 1. We shall now prove the following proposition that shows decomposability of A(x) in block 
diagonal form under the genericity assumption in Dcf. 2.1. 
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Proposition 3.2 Suppose A r in the series (3-1) is in Jordan canonical form 

If the Lidskii pscudovalues of A r _\ are generic in the sense of Def. 2.1, then there exists a formal power 
series 



Y(x) =Y + Y, x~ j Yj ■ det(Fo) ¥= , 

i=i 

such that 

B{x) = y- 1 ay -y- l y' , 

is a formal power series in block diagonal form according to the block decomposition of A r . 
The above is a consequence of the following 

Proposition 3.3 If the Lidskii pseudovalues of A r -\ are generic, then we can find a gauge H linear in 
x^ 1 such that the first two leading terms of 



have the following form 



A(x) = H~ 1 AH-H~' l H' , 



(0 ni )"-(0 ns ) = A r 



.4 



7—1,2 



\ 





V o 



A 



(3-12) 

(3-13) 
(3-14) 



r — l,s 



where A r —\ t i are square matrices of dimensions rii x n,; such that the only non-zero entries of A, — x,i « r e 
the yA r _ x,ij , that is, the bottom left hand corner. 

Proof. Let H = Hq+Hix -1 . We first partition Hq and Hi into s x s blocks according to the decomposition 
of A r . We will denote these blocks by H^} and respectively. The block is a rectangular matrix 
of dimension nj x n^. We shall similarly partition A r ^\ into s x s blocks and call these Ajk for simplicity 
(i.e. we suppress the index r — 1). 

The coefficients of x r and x r ~ 1 in (3-12) can now be written in the form 

A r Ho — H()A r = 

A r H\ — H\A r = A r —iHo — UqAf—i 

We can write these equations in block form, in which the jk th block of both sides are given by 



(0) o-(O) A/ - 



AfjH 



(3-15) 



i=i 



11 



where Afj is the rij x rij dimensional shift matrix 
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We want to find Hq and Hi such that (3-15) is satisfied with A r _\ given by (3-13). Let us denote the 
following linear operator from the space of Uj x nk rectangular matrix to the same space by Cjk 



Cjk ■ Mat(n.j, 7ifc) — > Mat(nj, n&) 



X^AfjX- XN k 



(3-16) 
(3-17) 



The first equation in (3-15) means that is in the kernel of Cjk, while right hand side of the second 
equation in (3-15) has to be in the image of Cjk- To prove the proposition, we will need to know the 
structure of the kernel and image of Cjk and show that under the genericity assumption in Def. 2.1, 
(3-15) is solvable. 

The following can be verified by straightforward calculation 

Lemma 3.1 The kernel of Cjk is of dimension min(rij,nk) and it is spanned by nj x n k matrices J of 
the form 



J 



/0 




Vo ... 

/Jl ... 

Ji 



J x ... 

Jj 







j 



Vo 



Jn k -1 
Jl 

/ 



J?lj-1 



Jl J 



n k < rij 



rij < n k 



(3-18) 



(3-19) 



where Ji are arbitrary constants. That is, the kernel is spanned by nj x nk matrices with top right hand 
corner of the form 



Ji ... J, 
Ji 



m — 1 J T7i 
• • ■ Jm—1 



\ o 







Jl / 



and all other entries zero, where m = min(nj,nk) 

The next lemma (left as exercise) characterizes the image of the operator Cjk- 

Lemma 3.2 The image of Cjk is of codimension min{nj, rife) and it is given by nj x nk matrices Q such 
that 



■p-\-i,i 



o, P = i,.. 



(3-20) 



where m = min(rij,nk)- 
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This means that the "dual diagonals" to the ones appearing in the characterization of the kernel of Cjk 
are "traceless". 

We now look at the term V jk = J2t=i AjiH^ in (3-15). To simplify the notation, we shall denote 
the entry in the bottom left corner of Aji by a(j, I) and we write Hy* in the form (3-18) denoting the 
elements on the diagonals by b? k and h? k = if i > min(iij,nk) or i < 0. (i.e., the Ji in (3-18)). Let 
rrijk = min(rij, rik) and consider the terms 



E^ 

i=i 



jk 

7i j —p-\-i,i 1 



P = i,. 



in V ]k , as these terms determine whether the right hand side of (3-15) is in the image of Cjk- 
Let Cjk = min(rij — n k , 0), a simple calculation then shows that 



, rrijk 



(3-21) 



i=l 



1=1 



where R k is a term that involves only h q for I = 1, . . . , s and q < p + Cj k and /i^ fe = if g < 0. Recall 
that Zi is the number of blocks with the same size n±, I2 is the number of blocks with the same size ni 2 , 
etc. Let H k be the column vector with entries 



: ) . =hf, i 



1. 



,n k , g = l, 



where Wi = Ylq=i lq an d Pi ^ s the biggest number such that n Wi > i. 
More explicitly H k is the column vector 



H k = 




(3-22) 



Now let y fe be the following column vector 



rqk 



(v k ) = s Tv q 



i = L 



,n fe , g = l, 
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That is, V k is obtained from the vector H k by replacing h qk by the corresponding 'sums of diagonal' in 
V jk . 

( Vntl \ 



ysk 




(3-23) 



If we denote the bottom left hand diagonal by the 'first diagonal', the next one by the 'second diagonal' 
and so forth, and call the sum of their elements the 'trace of the diagonal', then the i th block in V 
consists of the traces of the i th diagonals. 

We can now write the equations (3-21) for j = 1, . . . , s in matrix form 



U k H k = V k 



(3-24) 



From (3-21), we see that the matrix U k is of the form 







u" 



A- 
1 

* c/, fe 



u k J 



where ★ denotes some non-zero entries. The matrices U£ are given by 





o(l,l) 


a(l,w nk -{) 





u k = 


a(w nt -i,l) ■ 





















that is, U k is the wf_ l matrix that has the principal block submatrix of the Lidskii matrix of size uin k -i 
on its top left hand corner and zero elsewhere. 

We now consider the term A r ^i_k in (3-15). Let Xi be the bottom left hand entry of A r -i t i, that 
is, its only non-zero entry. We can now write this term in matrix multiplication form as in (3-24) 

U k H k = V k (3-25) 

where U k is the following matrix 

/ XJ\ 

U k = 



v 
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where U k is the diagonal matrix given by 

(uf) = 0,j<w nk+1 

v ' 33 

ipi)^ = Afc,j>W„ fc+ i 

The right hand side of (3-15) can then be written in the matrix multiplication form 

(u k - U k ^j H k = V k , k=l,...,s 

By lemma 3-20, the condition that the right hand side of (3-15) is in the image of the operator Cjk means 
that the determinant of U k — U k has to be zero. This determinant is given by 

n k 

det (u k - L> fe ) = Yl dot (lJ k - ufj 

i=l 

The first non-zero entry of U k is on the w nk ~i + 1 row while the first zero entry of U k is on the w Uk row, 
therefore, U k — U k is given by 





o(ll) 


a(l,w nk + l) 


a(l,w„ k -i) 





i - u k = 


a{wn k +1,1) ... 


a(w 7lk + 1, W nk + 1) - Afc . 


a(w„ k + l,w„ k -i) 







«K t -i,i) 


a(w nk -i,w nk + 1) 


■ a(w nk -i, W nk -i) - \ k 
















-A fc 7 



where / is the identity matrix of dimension (raj_i — m„ t -i) . 

Wc sec that the determinants det (u k — U k ^j arc just products of Afc and the pseudo-characteristic 
polynomials of the Lidskii submatrices. Therefore the genericity assumption in Def. 2.1 would imply 
that the determinant of U k — U k is zero and that the right hand side of (3-15) lies in the image of Cjk- 
Q.E.D. 

We can now prove Prop. 3. 2 as a corollary of Prop. 3. 3. We can assume that the genericity condition 
holds and that A r and A r _\ are given by (3-13) in proposition 3.3. 
Let Y be a power series such that 

oc 

Y(x) =Y + Y, x~ j Yj , det(Fo) + , 

3=1 

and 

B{x) = Y~ X AY - Y~ X Y' , (3-26) 
for some Laurent series B in x. The coefficient of x r ~ l in (3-26) then gives 

i 

A r Yi - YiA r = A r _{Yi-x - A r - X + ^ A r^~i - Yi-iB r -i - (r - i)li-r-i 

1=2 

where we have set for convenience Yk = if k < 0. 
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Assuming that B is in block diagonal form, we will show that (3-26) is solvable. We first write these 
equations into block form, as in the proof of proposition 3.3, in which the (j, k) blocks of both sides are 
given by 

~ y}k^k = A r ^Y^ 1] - Y^ 1] A r _ hk + ( ]T A r _iYi_i - (r - i)Yi_ r _A 

\l=2 / jk 

E^^r (3-27) 

1=2 

where Xjk denotes the jk th block of X. Provided that the right hand side is in the image of Cjk, (3-27) 
determines Y^} up to the addition of an element in the kernel of Cjk ■ We now are going to show that this 

arbitrariness in the solution of the i^ 1 equation can be used to guarantee the solvability of the (i + l) th 
equation, namely in such a way that the right hand side of (3-27) lies in the image of Cjk in the (i + l) th 
equation. 

Let J be a matrix in the kernel of Cjk as in (3-18), and A r ^i,j be the rij x rij matrix 

/ ... 0\ 
... 

\\j ... 0/ 

where Xj ^ 0. When the genericity condition is satisfied, we also have Xj ^ Xk if rij = rik, j ^ k. 

Recall that the bottom left hand diagonals are called the "first diagonal" , the next one 'second 
diagonal' and so forth, and the sums of their elements are called the "traces of the diagonal" . A simple 
calculation shows that the trace of the I th diagonal of Ar-ijY^. ^ — Y^ k 1 " ) A r _i J t is 

^jJ c {k.j)+l " ^kJc(j.k)+l 

where c(l,m) = min(ni — n m ,0) and Jk = if k < 0. Note that if the blocks have the same size then 
c(l,m) = 0, but then Xj ^ Xk due to our genericity assumption: if they have different sizes then either 
c(l,m) or c(m, Z) is strictly negative, and hence the system for the Jfc's is triangular with determinant a 
power of Xj (or Xk, depending on which one is the "long side"). In any case we see that the traces of 
the first rrijk diagonals can be chosen to be any value by a suitable choice of J, therefore one can always 
choose 1 such that the right hand side of (3-27) lies in the image of Cjk- Hence (3-26) is solvable 
with a block-diagonal B. Q.E.D. 

By application of Prop. 3.2 and Prop. 3.3 we can reduce to analyzing a single block of the final 
connection and we accomplish this final step in the following proposition. 

Proposition 3.4 Let A(x) be a n x n matrix of degree r — 1 and formal Laurent tail, such that the 
leading coefficient is Al plus a nilpotent matrix of rank n — 1. Suppose also that the Lidskii coefficient 
of the sublcading term (in the constant gauge in which A r is in Jordan canonical form) does not vanish. 

A{x) = x r - x A r + x r - 2 A r ^ + ... (3-28) 

Then there exists a nonsingular formal gauge Y(x) of the form 

Y(x) =C(l + 0{x- 1 )) , (3-29) 



A i — i,j — 
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such that the connection d x — A(x) is formally gauged to the connection d x — D(x) of the form 



n 


1 rn—1 

+ -E 

nx z — ' 

3=0 


t 3 W 


G 

X 


t rn := nX , G := diag(0, 1, 2, 


..,n 


-1) 




"00 





X 




1 








H = 


1 










o o '■ 













1 






(3-30) 



Proof. First of all we assume without loss of generality that A = since we can gauge it away by a scalar 
gauge transformation. Next, we can conjugate the connection by a constant gauge so that the leading 
coefficient becomes the nilpotent Jordan block TV. Then, as shown in Prop. 2.1, we can perform a second 
conjugation (constant in x) which recasts the problem in the form 



A(x) i-> A(x) = x 



r-l 



Ai 





A 2 
Ai 




A n - 















Ari-1 
A n -2 



A 2 
Ai 




* 
Ai 



(3-31) 



where the coefficients Xj are defined as the coefficients of the formal expansion of the eigenvalues of 
A(x) in Puiseux series of q — x x l n 



vM) = y(u J q) 



y(q) = q (r - 1)n J2 X ^~ j = ^ _1 EV 

3=1 3=1 



(3-32) 
(3-33) 



In fact we should identify the matrix x 1 r A(x) with the perturbation problem M(e) = M(l/x) used in 
Prop. 2.1. 

We use the same symbol A(x) for this new connection in "canonical" form to economize on notation. 
At this point we can perform the shearing transformation as done in Prop. 2.1: the difference is that now 
the transformation is a change of gauge and not merely a conjugation. Therefore we introduce another 
formal connection B which will be a formal Laurent series in q = x x l n (by choosing in an arbitrary but 
fixed way the determination of the root) 

B(q) := SAS- 1 + S'S- 1 
S :=q G , G:=diag(0,l,2,...,n 

Finally we change the variable of differentiation from x to q = x 1 ^ 71 



1). 



(3-34) 
(3-35) 



= nq n ~ 1 B(q)^ =: nB{q)^>. 



(3-36) 



where the matrix B(q), as a consequence of the shearing, has the following structure and enjoys the 
following properties. 



B(q) 



„-R-l 



X\C +-Br-\ 

q 



R 

E 



Brf 
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A n -1 



Ji-n 



1 ... 

10 

o '•. 

... 1 

1 

Bj = bjC 3 , bj diagonal matrices 
R = rn — 1 

(b R -j) u = ,j = 1, . . .n- 1, 1,. 

w5(w5) = siB^rr 1 , 



(3-37) 



SI := diag(l,LU,uj ,...,LO n ), a-> 



(3-38) 



Since the leading coefficient is nondegenerate we may apply the standard theory of asymptotic expansions 
as in [19, 18, 11] stating that we can find a formal solution in the form 

$ Z(q)c TM 

Z = Zo + Oiq- 1 ) 
R 

T(q) = ^-T^"+T ln(g) 



However it is more suitable to derive a slightly improved statement as in the following 
Lemma 3.3 There exists a (unique) formal solution satisfying the conditions 

R t 

T(q) = + t ln(q) , tj scalars 



Z(q) 
Z{oo) 



ST x Z{uq)C 

w , 



where W is the following eigenvector matrix of C 

cw = wn- 1 , W,, 



-(i-i)O-i) 



(3-39) 
(3-40) 

(3-41) 
(3-42) 

(3-43) 



Proof. 

This lemma is the counterpart of Prop. 2.1 in the setting of formal gauge equivalence. We start by 
remarking that the "periodicity" properties of Z and X" are compatible with the periodicity properties 
of B 



nB = nZT'Z- 1 + Z'Z' 1 
LuT'{ujq) =CT'{q)C- 1 . 



(3-44) 
(3-45) 
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The differential equation is conveniently rewritten as 



-Z' = BZ- ZT . (3-46) 
n 

Let us expand Z in (formal) Laurent series 

oo 

Z(q) = J^q-iZ j , Z Q := W . (3-47) 

3=0 

The periodicity of Z implies the following structure for the coefficients of the power series expansion 

Zj = & ZjW , Zj — diagonal matrices. (3-48) 
We first claim that any such series can be factorized as follows 

oo oo oo 

Y^(T j C j ZjW = ^2q- e C £ Pe W x^2q~ k j k n- k = (3-49) 

j=a e=o k=o 

oo / oo \ 

= J^q-'C'ptWxexpl J^Q^k^) (3-50) 
e=o \fe=o / 

Ti(pe) = , 70 = 1 , po = 1 , 7j = scalars . (3-51) 

Indeed, by comparing the two sides term-wise in the expansion one finds the recursion relations 

n 7j - Tr(^) , (3-52) 

3-1 

Pj = Zj -jjl-Y^Pklj-k ■ (3-53) 
fc=i 

Once the coefficients jj are known the coefficients S k in the exponential of last expression can be similarly 
defined in a recursive and unique fashion. 

If we plug this factorized expression in the equation and compare the terms of the same power we obtain 
for the coefficient of q K ~ 1 

K 

(B R _ 3 C K -ipK^W - tn^CK-ipK-jWn- 3 - 1 ) = , K = 0,...R. (3-54) 

3=0 

We recognize that this recurrence relation (for K < R) is exactly the same appearing in the proof of 
Prop. 2.1, which implies that t R = Ai. Since Br = t R C and C~ J W = WSV we obtain an equation for pk 

1 - 

Pk - C-'pkC = —J2 (tR-jC-i-W-jC^ 1 - C- K - 1 B R ^ J C K -^p K ^ J ) . (3-55) 
R 3=1 

This equation fixes the diagonal matrix px up to addition of a multiple of the identity (i.e. modulo the 
trace), but since pk is traceless, this equation suffices in fixing it unambiguously. Moreover this equation 
also determines t R -K from the fact that the LHS is by default a traceless matrix. Hence 



t R - K = ^Tr ^C- k - 1 B R -k + J2 {C^- 1 B R _ 3 C K -> p K _ 3 ) | . (3-56) 
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This also shows that tR-K = Xk+i since the above is the same recurrence relation that defines Xk, 
K < R. 

For the coefficients of q~ , K > 1 we have instead the following relations (we set t-j = for 
convenience in writing the following formula) 



R 



K 

n 



K-X 

C K p K W + - V (£ - K)S e ^ K C e p e C K - e W 

71 ' 



(3-57) 
(3-58) 



Multiplying both sides on the left by C K and on the right by W 1 we obtain 

K / 

tRPK+R - t R C 1 pk+rC = p K + 2J ( tj&PK+jC 3 - C K B j C k+:, pk+ ] 



j=-K 



K-l 



(3-59) 
(3-60) 



e=o 



This equation determines pk+ r. (modulo trace as discussed above) . Moreover tracing both sides gives an 
equation for 5k 

K8 K = -Tr ( C- K B J C K+ ip K+ )j , K > 1 . (3-61) 

This determines the scries Z uniquely Q.E.D. 

We can now conclude the proof of Prop. (3.4): to this end we gauge transform once more the 
connection d q — nB(q) by means of the formal gauge Y := Z(q)W~ 1 q G , namely 



d q — nB(q) i— > d q — nD(q) 

T'(q) 

G 



D(q) := q-°W {z'^BZ - ^Z^Z'j W 
The formal series D{q) enjoys the periodicity 

D(uq) = -D(q) 

OJ 

Restoring the independent variable to x we have the connection d x — D(x) where 

3 



(3-62) 
(3-63) 

(3-64) 



1 1 R 

D(x) = -qD(q) = —J2tj 



3=0 















X 


1 

















1 


































1 






R 



G 1 v — _ t , . i G 

— — > tjHixy — 

nx nx * — ' nx 

3=0 



(3-65) 



At this point we can summarize the chain of gauge transformations as 
Y~ X AY - Y- 1 Y' = D(x) 



(3-66) 
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Y(x) = q- G Z{q)W- 1 q G = Y j q- ] q- G C j z j q C 



3=0 



3=0 



X)Zj 



(3-67) 
(3-68) 



Since zq — 1 and zj are diagonal matrices, a direct inspection shows that the leading coefficient of Y(x) 
is 



n-l 



Y(x)= \l + J2^ jz 3+ <D ( x ~ 1 ) 

3=1 



(3-69) 



Therefore the formal gauge V(x) is nonsingular (note also that it contains only integer powers of x, in 
spite of the intermediate steps). Finally we claim that the constant term in Y(x) is in the ccntralizcr 
of A/" 4 : indeed by inspecting the leading term of D(x) we see that it is the same as the leading term of 
A(x) (in the "canonical" form 3-31), and hence the constant term of Y(x) must commute with it and 
so commutes also with its canonical form, which is Af. We can incorporate 1 + ^TV- 7 zj it in the first 
constant gauge, thus completing the proof. Q. E. D. Prop. 3.4. 

Remark 3.2 For later reference we point out that the upper triangular matrix U := Y^jZi^Zj is uniquely 
fixed once the connection is in the canonical form (3-31). Moreover, from the recursion relations defining yj's it 
is obvious that the entries of U are rational in Ai = t rn -i and polynomials in the other entries of the coefficients 
of A(x). 

Before proceeding with the analysis of isomonodromic deformations we point out the almost obvious 
Proposition 3.5 [Bare Isomonodromic Deformation] Introducing the matrices 



rn 

D°°(x) = — 5>W(fl 



G_ 

nx 



G := diag(0, 1, . . . , n - 1) 



T 7 C 



\x) 



3=0 

-jH(x) J , J = l, 



, rn . 



(where Ti(x) is the matrix appearing in eq. 3-30) we have the zero-curvature equations 



d 
dTj 
d 



J 



sj~oo,bare 



d 
dt K 



T, 



K 



d x -D°°(x) 



The kernel of all these connection is spanned by 



bare 



tg-o 



\J=1 



Proof. The above connection is nothing but the connection 

1 



q J V J , J=l, 



,7=0 



J 



(3-70) 
(3-71) 

(3-72) 
(3-73) 

(3-74) 
(3-75) 



This does not follow from the above expression directly because the diagonal matrices Zj are not multiples of the 
identity. 
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-the compatibility of which is absolutely trivial- written in a different gauge (and variable), namely by 
gauging with q~ G W. The kernel of these connections is trivially 



<fr(q) = exp 



j 



t In q 



Therefore we can set (multiplying on the right by the constant invertible matrix W 1 ) 

q-CWQW- 1 = q~ 



tybare _ Q -G w $ yy-l _ q -Gyy <j, yy-l q G q -G 



Since H{x) = q^WqClW^q , we find 



* bare (x) = exp 



^ 3 



Q.E.D. 



(3-76) 



(3-77) 



(3-78) 



Note that in Prop. (3.5) the sum extends to rn; namely t rn is nX, where A is the unique eigenvalue 
of the leading coefficient matrix. Among the times, each tj, j = Omod(n) is "trivial" in the sense that it 
could be gauged away by a scalar gauge transformation. However we prefer to keep them in view of the 
general case, where they will not be "trivial" anymore. 

In these considerations we have always placed the pole at x = oo: if the pole were at a finite point 
x = c wc should modify some of the formulas above in a trivial manner which is left to the reader to 
verify. The analog of Prop. 3.4 localized at a finite point is (we are not considering Fuchsian singularities 
in this section) 

Proposition 3.6 [Prop. 3.4 for finite poles] Let A(x) be a n x n matrix with a pole at x = 7 of degree 
r + 1 > 2 and formal Taylor series, such that the leading coefficient is Al plus a nilpotcnt matrix of rank 
n — 1. Suppose also that the Lidskii coefficient of the subleading term (in the constant gauge in which 
A r is in Jordan canonical form) does not vanish. 

A r A r -i 



^ ^ (x — 7) r+1 ' (x — j) 7 " 
Then there exists a nonsingular formal gauge Y(x) of the form 

Y(x)=C(l + 0((x- 7 ))) , 
such that the connection d x — A(x) is formally gauged to the connection d x — D(x) of the form 



tj 



rn 
1 1- , 



H J ((x- 7 )-') 



G 



Httx-j)- 1 ) 



nX , G := diag(0, 1,2,. 

" (x 
10 
10 



n(x — 7) 
..,n-l) 

7)- 1 












(3-79) 
(3-80) 

(3-81) 
(3-82) 

(3-83) 



Proposition 3.7 [Bare I-Defs for finite poles] Introducing the matrices 

tj G 



D\x) 



rn 



-H'ax-j)- 1 ) 



7^0 < x - 7) 

T]' hare {x) := -h"tt~ 



n(x — 7) 

~lY l ) J , j = l,...,rn 

C^' bare (x) := -D~<(x) , 



(3-84) 

(3-85) 
(3-86) 
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we have the zero-curvature equations 



" d 




dTj 


' d 


.7-7, bare 

1 J 1 U X 


' d 


_ C -f,bare Qx _ 1 


" d 


B 

'dt K ~ 





/y-7 ,bare 
1 K 



= 



2 K 



= 



The kernel of all these connection is spanned by 



exp 



(3-87) 
(3-88) 
(3-89) 
(3-90) 

(3-91) 



We conclude with the remark that the number of isomonodromic times in this situation is exactly the 
same as in the nonresonant case, and in fact this still holds for more general Jordan block decompositions: 
in a sense we may consider the cases under analysis "minimally" resonant. 



3.1 Isomonodromic Deformations 

The general situation in which we would like to define a set of isomonodromic deformations is that of an 
arbitrary rational connection A(x) with leading coefficients at each singularity consisting of (a conjugacy 
class of) an arbitrary Jordan form, under the suitable genericity assumptions on the Lidskii matrix for 
the subleading coefficient. 

At this point of our discussion the main difficulty is rather notational than conceptual: in view of 
Prop. 3.3 and Prop. 3.4 it should be clear that the approach can be "modular" by reasoning on each 
block in the local formal gauge at each of the singularities of the connection. 

It is probably an instructive warm-up exercise to consider the following simplified but nontrivial 
situation: let d x — A{x) be a polynomial connection of degree r — 1 with leading coefficient consisting 
of a single Jordan block of size n. There is no loss of generality in assuming that the connection is in the 
canonical form (3-31): this is tantamount requiring that the leading term is diagonal in the non-resonant 
case and fixes conveniently the constant gauge arbitrariness. More specifically we know from the proof of 
Prop. 3.4 that the coefficients appearing in the canonical form (3-31) are in fact the first n times, namely 
we will have 



A(x) = -x r 
n 



trn—n+1 
trn—n+2 



"rn — 1 



(3-92) 



where the parameters t rn -j, j = 0, ...,n — 1 are the first n times in the bare form (Prop. 3.5) of 
the connection. This can always be achieved -if necessary- by first putting the leading term in Jordan 
canonical form, then determining the parameters t rn -j and then applying a second appropriate constant 
conjugation. 
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Theorem 3.2 Let A{x) be a polynomial of degree r — 1 with leading coefficient in the gauge-fixed form 
(3-92). There exists a unique formal solution = A%> of the form 



=Y{x) 

rn 



* = (1 + U + 0(x~ 1 )) exp 



j 

U=l 



a;' "™ , (3-93) 



where U denotes a strictly upper triangular matrix uniquely determined, rational in t rn -i and polynomial 
in the remaining coefficients of Aj 's, independent of x and in the centralizer of TV. 
We then define 

Tj{x) := j (YH J (x)Y-') + , (3-94) 
where ()+ means the projection onto nonnegative powers of x. With these definitions, the operators 

8 X - A{x), dj - Tj(x) , J = 1, . . . rn (3-95) 
satisfy pairwise zero-curvature conditions, which ensure the Frobenius integrability of the Pfaffian system 

— * = A^> (3-96) 
ax 

d tj V = Tj^l . (3-97) 

The matrix U appearing in the formal gauge is a constant of all the motions and hence it can be disposed 
of by an appropriate conjugation by 1 + U . 

Note that although U is constant and in the centralizer of the leading coefficient it is not true that it 
can be set to zero a priori because setting it to zero is a (although conserved) algebraic constraint on the 
coefficient of the connection. We should rather think of U as "Casimirs" of the flows. 
Proof. The property of U being uniquely determined follows from Remark (3.2). We rewrite the Pfaffian 
system for Y as follows 

d x Y = AY-YD (3-98) 
d tj Y = TjY -Y Tj ,are , (3-99) 

where D(x) and Tj are are as in Prop. 3.5 for 

V bare := exp (j2 jK J ) ^ (3-100) 

is the solution of the Pfaffian system of Prop. 3.5. The definition of Tj implies that 

dtjYY- 1 = 0{x~ l ) , (3-101) 

namely the upper-triangular matrix U in the constant term is actually a constant of all the flows; note 
that the independent entries of U are n—1 since the matrix commutes with Af and depend polynomially 
on the (coefficients in x of the) entries of A, rationally on t rrl _i. 
We now want to verify the zero-curvature conditions; to this end one computes 



[d tj ,d tK ]Y ■ Y- 1 = djTjc - d K Tj + [T K , Tj] - [djT*? re - d K T h j are + [T* are ,T* are }) (3-102) 
[d tj ,d x ]Y -Y- 1 = djA-d x Tj + [A,Tj]-(djD-d x Tj are + [D,T} are }) . (3-103) 
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In the above equations the LHS is 0(x x ) since U is a constant of the motions, while the RHS is 
polynomial. Therefore both sides vanish identically. Q.E.D. 

Note that we could conjugate the connection Aix) by 1 + U since a posteriori this is a constant gauge 
of all the flows. However we could not gauge it away to begin with, since we did not know a priori 
whether it was a first integral of the isomonodromic deformations. 



Remark. We would like to comment about U being a constant in the proof of theorem 3.2. 

The situation is here very similar to the nonrcsonant case [11]: indeed the deformation equations that 
we have proposed here are not the most general that would preserve the leading coefficient at oo of the 
connection A{x). The residual freedom in this case as well as in the nonresonant case is by gauge action 
of constant (in x) transformations in the centralizer of A r . 

Suppose indeed we have a general solution {A, T} to the zero curvature condition (3-96) such that T 
is a polynomial. Let their formal solution be given by 



*= [1 + U + 0(x- L ) exp 



j=i 



(3-104) 



Then U need only be in the centralizer of A r (i.e. of M) but it could otherwise depend on the "times" 
in an arbitrary analytic way. Any such ^ can be obtained from multiplying a formal solution (3-93) by 
C = (I + U)(I + U)~ l . This has the effect of gauge transforming the solution {A°,T } by the gauge C, 
where {A°,T } is the solution in which U is fixed. Therefore any solution to the zero curvature condition 
can be obtained by a gauge transformation and there is no lost of generality in assuming U is a constant. 



3.1.1 Tau function 

We now turn our attention to a suitable definition of isomonodromic tau function. To this end we 
employ the strategy used in [6], namely of expressing the tau function in terms of spectral invariants of 
the connection A. This is motivated by the fact proven in loc.cit. that the standard Miwa-Jimbo-Ueno 
definition can in fact be expressed in terms of spectral invariants and this formulation is more suitable 
to be generalized to this case. 

First of all we make a few simple observations about the formal gauge Y used in eq. (3-93). Since we 
have 

A(x) = Y'Y- 1 +Y ( — YtjH{x) J - — W" 1 = 
\ nx nx I 

\ .7=0 / 

= Yx-iW ( ±Yt jX i&- X%W ~ lGWX ^ \ W-lggy-l + y/y-l , (3-105) 
\ nx nx I 

\ .7=0 / 

=: D(x) 

it follows immediately that the matrix P{q) := Y x~^W is an eigenvector matrix mod q-' rn - 1 + n where 
q = x~ . It also follows that the expansion of the eigenvalues yj of A coincides with the expansion of the 
eigenvalues of the bare system D{x). This leads to 

JTj = (YH J Y- 1 ) + = (xnYx-%WCL J W- 1 x%^ = (q J Pn J P~ 1 ) + + constant , (3-106) 

where the plus in the subscript always denotes the polynomial part in x. The constant appearing in the 
RHS is actually present only for J > rn — n + 1, namely for the deformations along the highest n times; 
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in any case the specific form of this constant (w.r.t. x) matrix is irrelevant for our purposes. 

Note that the expression in the RHS is nothing but a combination of the spectral projectors of the matrix 

A(x); if we denote by E a the diagonal elementary matrix we have 

n n 

pn-'p- 1 = J2 u^^PEvP- 1 = J2 u J{a - 1} n a , (3-107) 

(7=1 cr = l 

where II a is the rank-one projector on the eigenspace of the eigenvalue y a {q) (part of the cyclically 
permuted n-tuplet). By the properties of the spectrum of A we have 



U*(uq) = IL a+1 {q). (3-108) 



This ensures that the expressions 



q JJ2uj J((7 -^n a (q) (3-109) 

<T=1 

contain in fact only integer powers of x = q n . We can write the spectral projectors by the classical 
formula 

U a (q) = 1-—A^ a , (3-110) 

Tr(A - y a ) 

where the tilde denotes the classical adjoint (the matrix of cofactors). 

Note that once we have fixed the determination of the root q = and of t rn -i = (A,-) 1 /", there is 
a unique eigenvalue which admits the (in fact convergent) Puiseux series expansion 

VM = - £ *V- n + — + - E KH K q- n ~ K , (3-111) 
n y-^ nx n ^— ^ 

where this formula has defined symbols for the coefficients of the negative powers of the expansion. From 
an algebro-geometric point of view the cyclicity of the eigenvalues near x = oo means that the spectral 
curve has a branch-point of order n in P 1 x P 1 at y = oo even after desingularization; put otherwise the 
local parameter at near x — oo is q = x* . This means that we have 

tj = - res (y^x- J/n u) Ja y a+1 (q)\ dx = resx~™ydx , 1 < J < rn (3-112) 
i = - res I > ' y a (q) ) dx = - res Tr{A)dx =t - =t - ^— - (3-113) 




Hj = res x J/n LU Ja y a {q) | dx = ^-lesx^ydx . (3-114) 

./ rr — nci \ < ^ I J 



a=Q 



We remark that the first residue- formulas for tj and Hj are taken on the x-plane (the "base-curve") 
and they make sense because the expressions in the brackets contain only integer powers of x due to 
the cyclicity of the eigenvalues: the second residue-formulas are taken on the spectral curve around the 
point Coo (with local parameter x~ x / n ) which projects down to x = oo. Since the positive orientation of 
a small circle around Coo is the opposite than the positive (counterclockwise) orientation in the x-planc, 
this explains the difference in sign. Note also that the coefficient to of the bare system is not the residue 
of the eigenvalue but it is shifted due to the fact that the term G/x contributes to the residue (3-113). 
We are now ready to define the tau function for this guide-example 
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Proposition 3.8 The following differential is closed 

rn 

dlnr = ^2 H K dt K . (3-115) 

K=l 

Before entering in the (easy) proof let us comment on the motivation and shape of the formula. The 
reason for the name of this differential as "isomonodromic tau function" resides in the fact that for the 
nonrcsonant case [6] a similar formula involving residues of the spectral differential ydx on the spectral 
curve of the connection is shown to coincide with the definition in [11], which is of different nature and 
-on the face of it- not a spectral invariant of the connection. Secondly, this sort of expressions appear 
throughout the literature of dispcrsionlcss integrablc hierarchies, Seiberg-Witten models etc., we cannot 
find in the literature any statement or implication that the Miwa-Jimbo-Ueno tau function is defined by 
or equivalent to this ubiquitous formula. 

On a more technical point we also remark that the most similar setting is that of the so-called universal 
Whitham hierarchy [13], where the data arc an algebraic curve and a meromorphic differential plus some 
decorative data of local parameters near punctures. 

While the main characters appear similar (we have a spectral curve, we have punctures and local 
parameters and a meromorphic differential ydx) there are some significant distinctions regarding the 
coordinates on the phase space. In fact in the Whitham setting the periods of the differential ydx would 
be treated as coordinates on the phase space independent of the other "times" . In this case the periods 
of the differential ydx are not independent and in fact mainly uncontrollable. What is fixed instead 
are the parameters of formal monodromy and the Stokes' parameters, which cannot be recovered by 
inspection of the spectral curve alone. 
Proof of Prop. 3.8. 
We compute the closure of the differential 



dt,H K = ^vesx K / n d tj ydx = ~ T ( T\ res x f ^ K{a ~ 1] d tj Va{x) J 



dx . (3-116) 



In order to compute the variation of the eigenvalue we recall the classical formula describing the variation 
of a simple eigenvalue y of a matrix A under an infinitesimal deformation A i— > A + 5 A 

Sy = }- Tr (a^v 5a) . (3-117) 

Tr(A-y) V > 

In our case we should use the formula with S A = Tj + [Tj, A] because on the path of integration the 
eigenvalue is "uniformly" simple. We thus obtain 

- d tj H K = I res (XX^Ar ( ^ZJ^ x) (d x Tj + [Tj,A]))) dx = (3-118) 



Kb=l 



s Tv(A-y b )(x) 



= ™» E^ 6 - 1 ^*" Tr ( n ^)£ (E^^IW ) ] d* = (3-119) 

I ( n — - 

= res res V ^(b-D+Jia-i) x " ' Tr (n 6 ( g )n„(p)) ) dzdx . (3-120) 

JK x=ooz=oo I * — ' (X — Z) 1 

At this point the formula is almost obviously symmetric: attention is to be paid as to whether the order 
of residues can be interchanged. To see this, consider the kernel 

n ba {x,z) := -^—Ti{U b (x)U a (z))dxdz . (3-121) 

(x — zy 
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To conclude that the order of residues can be interchanged it is sufficient to show that there is no residue 
on the diagonal x = z. To this end we note that 

Tr(n Q (a;)n b (2)) = S ab + 0((x - zf). (3-122) 

Indeed the fact that Tr(II a (z)II&(x)) — 8 a b is nothing but a statement of simplicity of eigenvalues (x ^ oo): 
by continuity it holds also at x = oo. Next, we note that if C(x) is a locally diffcrcntiable matrix of 
eigenvectors, we have 

Tr(Ti a (x)Tl b (x)) + 0{x - z) = 
Ty(n a [C'C-\n b ])+0{x-z) = 

Tr(C , 'C , - 1 [n a ,IL J ]) + 0(x -z) = 0(x - z) , (3-123) 

where we have used the cyclicity of trace. This proves that the kernel tt a b(x, z) has only a double pole 
without residue on the diagonal (and only for a = b) and this is sufficient to have independence of the 
order of the residues. Q.E.D. 

Note that in fact the kernel introduced in the proof has some interesting properties because it is 
analytic on the whole spectral curve except on the diagonal where it has a rcsidueless normalized double 
pole. This is almost the same as the fundamental bidifferential of the spectral curve, the only shortcoming 
being the absence of a definite normalization around an isotropic basis in the homology of the curve. 

Before taking on the case of a general rational connection, we highlight the overall logic used to 
construct the Pfaffian system, prove Frobcnius intcgrability and construct the tau function: 

1. we start with the bare zero-curvature equations in Prop. 3.5 

2. we have dressed by the formal nonsingular gauge Y 

3. we have defined the tau function in terms of residues of a canonical differential on the spectral 
curve. 

The key of the proof of zero-curvature is the usual argument that the curvature of the Pfaffian system 
is on one hand a priori a polynomial expression in x whereas on the other side it is also a Laurent 
series, hence concluding that it must identically vanish: this argument -of course- hinges on the fact that 
the formal gauge Y is nonsingular. The deviation from the usual situation is not in the way we prove 
the compatibility, but rather in the precise form of the bare system of equations which was not known 
beforehand and which in the nonresonant case is in diagonal form. 

4 Fuchsian resonant case 

Isomonodromic deformations of resonant Fuchsian system in full generality have been addressed in [7] so 
that here we mainly collect some known facts for the reader's convenience. Consider a formal Fuchsian 
singularity (at z = 0) 

n*) = !+£ F ^ ■ (4-1) 

j>0 

We assume without loss of generality that A is in Jordan canonical form. Since the notion of resonance 
for (formal) Fuchsian singularities is that two eigenvalues differ by a nonzero integer, we split A into 
blocks of (Jordan blocks of) eigenvalues differing only by integers (a bouquet). It is not difficult to show 



Ti-(H a (x)U b (z)) 

x — z 



28 



that there is a formally analytic gauge in which the connection is completely split into block diagonal 
form according to the decomposition into bouquets. Hence we can assume that A has only eigenvalues 
differing by integers (i.e. consists of only one bouquet of eigenvalues). By a scalar gauge transformation 
we can actually shift the eigenvalue to zero. 

At the end of these simplifications we have a connection in which A is in Jordan canonical form with 
only integer eigenvalues which we assume in decreasing order 



n > m > n 2 > 



> n K 



(4-2) 



Each eigenvalue has a certain algebraic multiplicity and a certain geometrical multiplicity (i.e. the rank 
of the eigenspace), which is of no particular interest to us. 

We now consider a formal analytic gauge Y(z) = J2j>o^j z '' ana - see ^ ^ ne mos t "canonical" form 
under formal gauge equivalence 



Y' = FY - YB 



* = 7 + E*i 

j>a 



Y=l + 0(z) 



(4-3) 



Our goal is to have B as simple as possible. Writing out the coefficients of the power z we have 



fe-i 



kY k = [A, Y k ] +F k -B k + J2 (FjY k -j - Y k ^Bj) 



(4-4) 



We rewrite this as 



fc-i 



(fcld - ad A )Y k = F k -B k + J2 (FjYk-j ~ Yk-jBj) 



(4-5) 



The linear operator C k := (kid — cicIa) on the space of matrices nxnis invertible provided that k is 
not in the spectrum of ad\, i.e. provided that no pair of eigenvalues of A differ by k. If C k is invertible 
we can impose B k = since the solvability of the recurrence relation in terms of Y k is guaranteed. If 
L k is not invertible then B k must be chosen so that the RHS of (4-5) is in the image of C k . It is not 
difficult to see that the image of C k consists of arbitrary matrices with a zero block in the (j, £) block 



such that 



k. Therefore B k can be chosen to be zero in the complement of that block and it is 



then uniquely determined by (4-5) itself. By finiteness of the number of eigenvalues of A we can assure 
that C k is invertible for k large enough, namely that only a finite number of B k may need to be chosen 
nonzero. At the end of this procedure we always obtain a connection in the form 



B 



" n + V 


■kZ n °~ ni 


■kZ n °~ ™ 2 


-kZ n °~ n3 




^ z n0— n K 





ml + v 




*z ni ~ n3 












n 2 l + V 


~kZ n2 ~ n3 


























n K -il + V 


^. z ni<-i—n K 












n K l + X7 



(4-6) 



where V denotes a nilpotcnt matrix in Jordan canonical form and * denote the only possibly nonzero 
coefficients and constant in z. Each diagonal block in the above decomposition has dimension equal to 
the algebraic multiplicity of the corresponding integer eigenvalue of A. The block diagonal shearing 



{z n °l,z ni l,z n n,...) 



(4-7) 
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(here each identity is of the appropriate dimension) recasts the connection to a nonrcsonant one with the 
only eigenvalue zero of the simple form 



-G 



Bz 



G 

z 



T 

z 



T := 



V 


★ 













V 
















V 




























V 














V 



(4-8) 



Note that T is constant. The *'s in (4-6) and (4-8) represent the same coefficients and are defined up to 
action of constant gauge transformations in block diagonal form, each nonzero block of which consists of 
the centralizers of the diagonal blocks of F. The orbit under this ccntralizer group is what defines the 
local monodromy and hence must be preserved by the isomonodromic deformation. 
The solution of last system is 



QOare _ z ± _ 

In other words a (formal) solution of the original system is 

* = Y{z)z 



G Z T 



(4-9) 



(4-10) 



where G is the diagonal matrix of integers used in the shearing (inducing the grading) and in general 
does not commute with T . Since T is upper (semi)triangular it can be put in Jordan canonical form T can 
by an upper triangular matrix P; 

* = Yz G Pz T — P- 1 (4-11) 

Since z G Pz~ G is analytically invertible, it can be reabsorbed in the definition of Y, so that without 
loss of generality we can always assume the formal solution (4-10) to have T in Jordan canonical form 
and hence we will denote T by J in the sequel. At this point the residual arbitrariness of this (formal) 
solution is by multiplication on the right by a z-independent matrix S in the centralizcr of J and such 
that z G Sz~ G is analytic; this implies that S must be at the same time in the centralizer of J and in 
block upper-triangular form according to the minimal decomposition of G into blocks which are multiple 
of the identity matrix. 

To put it in a different way, S must be in the intersection of the positive root spaces of Ad z G with 
the centralizer of J. 

Summarizing this discussion and restoring the generality of all the steps we have obtained the 

Proposition 4.1 [Formal solution for resonant Fuchsian singularities] Let A = ^f- + 0(1) be the matrix 
of a Fuchsian singularity at z = 0. Then there exists a (formal) solution $ of the form 



(4-12) 



Y(z)z G z J 
where 

1. Y(z) is analytically (formally) invertible at z = 

2. J is in Jordan canonical form with distinct eigenvalues with real part in the interval [0, 1) 

3. G is an integer valued diagonal matrix such that -within each block corresponding to the same 
eigenvalue of T— the integers form a weakly decreasing sequence which distinguishes the eigenvalues 
of the same bouquet. 



4. The spectrum of Aq coincides with the spectrum of G + J and moreover A$ = F (G + J)Yq 1 
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Such solution is unique up to ordering of the eigenvalues of J and by multiplication on the right by a 
^-independent matrix lying in the intersection of the centralizer of J and the nonnegativc root subspacc 
of Ad z c 5 . Moreover the gauge Y(z) is actually a convergent series if A(z) is convergent in a punctured 
disk around the regular singularity [19]. 

4.1 Isomonodromic deformation of resonant Fuchsian singularities 

Here we are rephrasing part of the content of [7] . Suppose we have an isomonodromic family of resonant 
Fuchsian connections with poles at points jj, j = 1,.. . and let $(257) be the fundamental solution of 
the family. At each pole jj and by virtue of the previous discussion culminated in Prop. 4.1 there is a 
z-independent nonsingular matrix Cj for which 

$ = Y 3 {z){z - 7j ft (z - yy'Cj . (4-13) 

The matrix Cj is defined modulo the group described in Prop. 4.1. The monodromy is then 

Mj = f , c 2 ' 7 '' C, . (4-14) 

We see that under a monodromy preserving deformation the matrices Cj can vary arbitrarily by left 
multiplication of a matrix in the centralizer of Jj (which is the same as the centralizer of e 2l ' Kji since the 
eigenvalues of Jj by construction do not differ by integers). 

Under a continuous monodromy preserving deformation we have near each pole jj 

= 7, Yj^jj-Yr 1 + Yj(z jjftCjCj-'iz - jj)- ^-' + YjYfi . (4-15) 

The last term is analytic at 7j . The first term has a simple pole and is the standard term in Schlesinger 
deformations. The second term may have poles of higher order: indeed CjCj needs only to belong to 
the centralizer Lie algebra of Jj but not necessarily to the centralizer of Cj nor to its nonnegative root 
subspace. A case in which the second term is certainly absent is when the centralizer of Jj is the Abclian 
algebra of diagonal matrices, which corresponds to the case of nonrcsonant Fuchsian singularities; another 
case is when Jj has only one Jordan block, for in that case the centralizer is upper triangular and hence 
certainly in the nonnegative root space of Ad z G . In case Jj contains more than one irreducible block with 
the same eigenvalue then the centralizer is not upper triangular (see Figure 1 for representations of such 
a centralizer) and hence conjugation by (z — r yj) Gj may have poles of higher degree (at most the index of 
resonance, i.e., the maximum integer difference between two eigenvalues). 

This situation is the most general as explained in [7]. We can regard the dependence of Cj on the 
deformations as "pure gauge" in the sense that it is arbitrarily defined and can always be disposed of by 
a rational gauge equivalence which does not move the position of the poles, or -which is the same- by 
solving a certain Riemann-Hilbert problem. 

Suppose indeed that $, are two isomonodromic families of a resonant Fuchsian connection with the 
same monodromy representation {Mj} and with same residual spectrum at the Fuchsian singularities 

Sp( res $' $~ x ) = Sp( res $' -I" 1 ) , Vj . (4-16) 

Z=~tj Z=-fj 

This implies that the matrices Jj (the Jordan form of the monodromics) and Cj are the same for the two 
families. This also implies that 

G(z)-.= §{z)$- x {z) (4-17) 

5 The operator Ad z o splits GL(n,C) into eigenspaces with eigenvalues z u ~ G ii ; the direct sum of the subspaces with 
nonnegative exponent for this eigenvalue is our definition of nonnegative root subspace. For example if the entries of G are 
strictly decreasing then the nonnegative root subspaces are the upper semitriangular matrices. 
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Figure 1: A nilpotent matrix in Jordan canonical form and the shape of its centrahzer: the oblique 
segments represent entries with the same numerical value [1]. 



is a matrix function defined on the punctured plane and single- valued. One sees from the asymptotic 
representation of $ and <J> that 

G{z) = %{z - jjfiCjCr^z - vr^Yj + 0{z - 7i ) (4-18) 

where we have used that Hj := CjCj is in the centrahzer group of J 3 because C, l e 2 ™ J ( ', = Mj = 

Mj = Cj~ 1 e 2l '* Jj Cj . This shows that G(z) has at worst poles at the jj's and hence it is rational. 
We can think of G(z) as the solution of the following RH problem 

G + {z) = G-{z)Y 3 (z~ 1 jfiH ] {z- 1] )- G iYf 1 , z e {\z - 7i | = e} (4-19) 

where G-{z) is analytically invertiblc in each disk and equal to Yj(z)Yj~ (z) (recall that both Yj and Yj 
are actual convergent series around jj [19]). Since the jump matrices arc analytic in each punctured small 
disk, the "exterior" part of the solution of this RH problem can be analytically extended to a rational 
function G(z) in the punctured plane CP 1 \ U{7j}. 

Vice- versa suppose that given an isomonodromic family $ we want to pass to another isomonodromic 
family $ with the same monodromy matrices Mj and same residual spectrum (4-16). To this end we 
should find a solution to the Riemann-Hilbert problem (4-19) with preassigned arbitrary matrices Hj 
depending on the deformation parameters (but independent of z) and in the centrahzer group of Jj. 
The "exterior" part G+ of the solution of the RH problem defines by analytical continuation a rational 
function G = G+ which transforms the family $ in a family $ with Cj = HjCj. The solvability of this 
RH problem can be assured by the argument which will be used later in the more general setting of Thm. 
5.3 (see Remark 5.1). 

This discussion means that for resonant Schlcsinger systems the Schlesinger equations are still con- 
sistent and any other isomonodromic family is obtained from a solution of Schlcsinger equations by a 
rational gauge equivalence constructed from a solution of (4-19). 

5 General Case 

We now address the full-fledged general case. As mentioned earlier the conceptual difficulty is an a par 
with the one involved in the Isomonodromic deformation system studied in Thm. 3.2 but we need to set 
up a good deal of notation. 

Let T> be an effective divisor on CP 1 

= ^00 + ^(^ + 1)7, roo^l, r 7 >0. (5-1) 

7£C 
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It is understood that only a finite number of points 7 appear in the above sum. We now consider an 
arbitrary rational connection of the form 



A:= y, ^» ( 5 - 2 ) 

AM:=E(^7TT (5-3) 



j=o 



A 0O (x):=Y l xJ ~ lA co,J ( 5 " 4 ) 
j=i 

where the leading terms of the singularity at each irregular singularity (r 7 > 0) may have an arbitrary 
Jordan canonical form provided that the Lidskii pseudo-eigenvalues of the second leading term subordi- 
nated to the Jordan form of the leading term do not vanish and are distinct. No restriction is assumed 
on the residues at the Fuchsian singularities. Note also that we are assuming that there is at least one 
irregular singularity and we have placed it at 00. 

5.1 Generalized monodromy data for resonant singularities 

We now proceed along the lines of [18] to define the generalized monodromy data for our system of linear 
ODE 

dtf> 

— = A^ (5-5) 
ax 

Let s 7 be the number of blocks in the Jordan form of ^4 7 ,r 7 and nj be the dimensions of the blocks 
for a non- Fuchsian singularity (r 7 > 1). Denote by f? 7 be the diagonal matrix 

t 1 t 1 - 1 t 1 - r? 7 4- 1 t 1 t 1 - 1 t 1 - T)T + 1 



r> i- I "0.1 "0,1 x "0.1 '"1 ~ x "0,s "O.s 

Rry = diag — ^ T 

\ n{ n{ 

where k are scalars. 

We can now use theorem 3.1 and the results in Section 4 to construct formal solutions of (5-5) near 
each pole. 

Proposition 5.1 Let A be given by (5-4) such that at each pole 7 € T>, -A 7jT . has the Jordan normal 
form 



[A,r T ] = (A£)-"(X&) 



and that the genericity condition in definition 2.1 is satisfied at each pole if r 7 7^ 0. Then for each pole 
such that r 7 ^ 0, there exists a unique formal series Y J 

00 

F>-7) = r 7 +^(.T- 7 )^ 7 , det(V)^0 (5-6) 
3=1 

such that 

V =r 7 exp(Q 7 (x-7))(a;-7)^ (5-7) 
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is a formal solution to the ODE (5-5), where Q-y(x — 7) is the matrix 



/ r 1 W j 



rn 2 njj 



"j>2 



7,s 



3,' j J 
3 = 1 J ' 



where 7i 7 .i are nj x nj matrices defined as in (3-30) and are scalars. 
When r 7 = 0, there exist a solution of the form 

^ = r 7 (.T- 7 ) G ^(x-7) J7 . 



(5-8) 



(5-9) 



where G 7 is a matrix determined by the integral differences between the eigenvalues of A^„ as in section 
4, Y 1 is a convergent series as in (5-6) and J 7 is a constant Jordan matrix. 

The uniqueness of Y 1 can be seen by first making a gauge transformation of A such that A^ >r is in 
Jordan canonical form. 

We can now use this formal solution to define the monodromic data of the ODE (5-5). We have the 
following 

Theorem 5.1 There exists a finite number of asymptotic sectors with vertex at x — 7, such that 
{5. 7 } * s a covering at x = 7. Moreover, there exist fundamental solutions of the linear ODE 

— = A$ 

dx 

asymptotic to the formal series solution of the form (5-7) or (5-9) near x = 7. 

The proof of this result can be found in [18] or [19]. Such solutions ^ 7 are analytic outside of the divisor 
V and on the universal cover of CP 1 \ V. 

In particular, we can now cover the punctured Riemann sphere CP 1 /!? with sectors Sj, 7 £ "D and 
define the stokes data by using the asymptotic form (5-7) and the sectors Sj. 

We have the following 

Definition 5.1 Let {S?} be a covering of the punctured Riemann sphere CP 1 /!? that satisfies the con- 
ditions in theorem 5.1 and let be the fundamental solution that is asymptotic to the formal solution 
of the form ( 5- 7) in S? . Let 



We shall call ^S^,C^f^ a Stokes phenomenon of the ODE (5-5). This, together with the set of variables 
tg j at each pole x = 7, r 7 7^ 0, and the J 1 for r 7 = 0, is called the monodromic data of the ODE (5-5). 

Here the matrices C£j contain the usual definition of Stokes' matrices and connection and monodromy 
on the same basis (and redundantly). 
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5.2 Isomonodromic deformations 



In view of Thm. 3.2, at each singularity 7 € I? we can find a formal nonsingular gauge Y 1 that gauges A(x) 
to the localized version of the bare form D^(x) advocated Prop. 3.6. The bare connection/deformation 
at x = 7 is now the direct sum of as many bare systems of the form given in Prop. 3.7 as the number of 
blocks in which the system has been decomposed. All these bare deformation/differential equations are 
then dressed by the same dressing matrix Y 1 . 

In order to be more explicit consider an irregular singularity x = 7. Suppose the Jordan form of the 
leading coefficient at x = 7 is as in Prop. 5.1. In symbolic notation we denote by t 7 the times "attached" 
to the pole at x = 7. The number of them is 

#{t 7 } = rn , (5-10) 

namely the same as if the system were nonresonant. 

The following theorem allows us to consider the direct sum of many "bare" Pfaffian systems of the 
form in Lemmas 3.5, 3.7 and obtain a solution of another Pfaffian systems after a "dressing" using a 
formal gauge equivalence. 

Theorem 5.2 Suppose A is a matrix-valued function given by (5-4) andY 1 is a formal series holomor- 
phic and invertible at 7 



y^(x - 7) = v + E(- T - riV > det ( r °) * 



i=i 

OO 

Y°°{x) = F °° +Y^x~ i Y9° 

such that 

B~<{x) = (Y^y 1 AY~< - (FT 1 (Y~<y , 

is a formal power series in (x — 7) in block diagonal form according to the block decomposition of Aj^ r , 
and similar expression for 00. 

Let T 7 (x) (T°°(x)) be matrix-valued 1- forms polynomial in (x — 7) _1 (x, respectively) such that 

[d x -B^{x),d-T\x)\ = 0(x-j) (5-11) 
[d x -B°°(x),d-T°°(x)] = Oix- 1 ) . (5-12) 



Then the Pfjaffian system 



[d x - A{x)] * = (5-13) 
[d-M(x)]V = (5-14) 



is Frobenius compatible, where M{x) is the matrix-valued rational 1-form such that the singular part of 
M. near each pole is given by 



M 



(x) = (yTfOz;- 7) (FT 1 ) +0(1), x^ 7 ^oo (5-15) 
where X pp denotes the principal part of X at x = 7. Near x = 00, M. is given by 

M(x)= (Y 00 T 00 (x)(Y oo y 1 +dY°° (Y 00 )^ + +0(x- 1 ), x -> 00 (5-16) 
where X + is the polynomial part of X and Yq 00 is some analytic function in the deformation parameters. 
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Proof. Let M.(x) = Y^, M. v {x)dt v : here the label v is a generic label used for the collection of all times at 
all singularities. According to this definition each M. v is singular only at one of the poles 7 G (D). Let 
v be the pole of M. v {x). A simple calculation shows that 

[d v ,d x ] Y v ((rr) _1 = duA - d x M v (x) + [A,M v [x)\ - {d v B v - d x T v v + [B v ,T»(x)}) 

where M v {x) = {y v T v {Y v )~ 

By using similar argument as in the proof of theorem 3.2, we see that 

[d x -A{x),d v - M v {x)] = 0{x-v), x^v 
\d x -A{x),d v -M v {x)] = 0{x- r ), x^oo 

Therefore, by Liouville's theorem, we have 

[d x -A(x),d u -M u (x)]=0 

Similarly, we have 

[d v - M v (x),d,j, - M^(x)] = <D(x-v), x^v = n 
[<9„ - M v (x), <9 M - M^x)] = 0(x~ 1 ), x—>i> = fi = oo 

when the pole of M. v (x) and M^x) are the same point. 

Now consider the case when the poles of M.„(x) and M.^{x) are not the same. Then near x = v 1 we 
have 

d v Y" = M v {x)Y v + Y U T» 
dpY v = M^V 

We also have similar equations with v replaced by (i. From this we have 

[d v , d M ] Y v (YT 1 = d v M v {x) - d^Mvix) + ^(1)^,(1)] 

near x = v and similar equation with [i replacing v. By using similar argument as before, we see that 

d v M v (x) - d^M u {x) + [M^{x),M v {x)\ = 

this concludes the proof of the theorem. Q.E.D. 

We can now consider isomonodromic deformations of the ODE 

— = A$> 

ax 

Let {S£, be a covering and fundamental solutions as in definition 5.1 and let |c£^, J 7 , fc | be the 
monodromic data associated to it. The isomonodromic problem is the following. Given {5^,4'^}, how 
should one deform the solutions such that the monodromic data |c^, J 1 fc | is fixed? 

The isomonodromic problem is only defined when a covering and solutions {S^^^} is chosen and 
should be thought of as deformations of the fundamental solutions. For if the monodromic data of 
{S^,^^.} is fixed under deformations, then by multiplying vf^ to the right by matrices CJ? that depend 
on the deformation parameters t, one obtains a pair {S£ ,^^C^} in which the monodromic data is no 
longer fixed. 
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We can now prove the following theorem, which classifies the monodromy preserving deformations 
when the genericity condition in definition 2.1 is satisfied. 

In the following we denote the matrix-differential one form of the bare deformations as 

(ml rnj^ \ 

c7d 7 + ]T i? , ' 6ore dt} ilJ . . . ,c7 7 d 7 + T ^; bare ^U ( 5 " 17 ) 
.7=1 J=l J 

where the notation is as in Props. 3.5, 3.7. 

Theorem 5.3 Let A be a rational matrix-valued function given by (5-4) and let T> be the divisor of poles 
such that r 7 = for 7 £ V and Di be the divisor of higher order poles. Suppose the genericity condition 
in definition 2.1 is satisfied at each pole 7 € T>\ and that 00 € T>\. 

Let {5^,^^} be a covering and fundamental solutions that satisfies the conditions in Def. 5.1. Then 
the monodromic data of the ODE (5-5) defined by , ^j?} are preserved if and only if satisfies the 
differential equations 

&n = (5-18) 

where T is the 1-form 

T = (y^T^'^iY 00 )- 1 + dY °° (T^r 1 ) + ^ (y"'T"'' bare (Y"')- 1 ) pp + 

+ y A 7 . (x - 7 )- x d 7 + (Y^{x - jf-dH^H^ix - 7 )~ G - {Y-ty 1 ) 

where T 7 ' are is defined in (5-17) and Y 1 , G 7> and Y 00 are given in proposition 5.1, while H 1 is a matrix 
in the centralizer of J 7 which is constant in x and depends analytically on the deformation parameters. 
We use X pp to denote the principal part of X around the corresponding pole and X+ to denote the 
polynomial part of X . The exterior derivative d in the above expressions denotes derivatives with respect 
to the t in proposition 5.1 and the position of the poles 7, but with dtg k = 0. The dependence of Y" °° 
and H 1 on the parameters, that is, dy o °° and dH 1 , can be arbitrary, provided that H 1 remains in the 
centralizer of J 1 . 

Proof. The proof is essentially the one in [11]. By using similar argument as in the proof of theorem 
5.2, we see that in the case of r 7 = 0, we have 

[c^a,] Y^ (Y-*)- 1 + Y 1 [d*,^] {(x - 7 ) G -£T) (V(x - 7 ) G -^) _1 = d^x) - a 7 A(x) + [A(x),T 1 (x)] 

near x = 7 where 7^ is the coefficient of d7 in T. The left hand side is a positive series in x — 7 since 
the second term is a well-defined function of x and the deformation parameters, and Y 1 is a positive 
series in x — 7. By considering the singular behavior of the right hand side near each pole and then 
apply Liouvillc's theorem, one sees that the right hand side must be zero. Similarly, one can show the 
commutativity between the parameters t, 7 when one of the parameters involved is the position of a 
simple pole. 

For the cases that do not involve derivatives of the position of a simple pole, we can apply proposition 
3.7 and theorem 5.2. We see that the Pffafian system 

f \d x - A(x)] * = 

I 1 (5-19) 
1 d-Ti f = 
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is integrablc. 

Let \E'£ and be solutions of (5-19) as in Def. 5.1, where a, f3 £ D and may not be the same point. 
Since ^ and "ff solve the same equations (5-19), we have 

dn (nr 1 = d*f = d*^ (nr 1 + n*cg (eg ) _1 o^r 1 

hence dC^j 3 (c^j 3 ) = and the monodromic data is preserved. 

Conversely, let be deformed in such a way that the monodromic data defined by {S^, 1 ^^} is 
preserved. We see that 

an o^r^d^f (*£ 

for a, (3 E T>, where a and /3 may not be the same point. Therefore d^ (^P^) -1 is a globally defined 
meromorphic 1-form 

M(x)~dnmr 1 ■ (5-20) 

Its asymptotic behavior near each pole are given by (note that dQ 1 = T' y > bare ) 

d&ZQ&Z)- 1 = d*f hi) 1 ~ dV 1 (Y^y 1 +Y^dQ 1 {Y^r 1 

= (r 7 dQ 7 (y 7 ) _1 ) + 0(1) x^ 7 , r 7 >0 (5-21) 

V / pp 

d^(^)- 1 = A 7 , (.T- 7 )- 1 d 7 +(^(.T- 7 )^d^(^)- 1 (.T- 7 )^(^)- 1 ) +0(1) 

x -> 7 , r 7 = (5-22) 

for all a, /3 and 7 <G 2?, where X ~ W means that X is asymptotic to W near the corresponding point. 
One can also write down a similar equation for x — > oo. 

The equations (5-21) determines the 1-form ./Vf up to the addition of a constant in a;. 

To determine this constant, we can look at the behavior of M. near x = oo. This is given by 

dr°° (y 00 )" 1 + y^dQoo (y 00 )- 1 



We see that Ai has the form of T in the theorem. Q.E.D. 

Remark 5.1 The essence of Thm. 5.3 is that the dependence of Y°° and H 1 can be regarded as gauge 
arbitrariness, and fixing it yields a consistent Pfaffian system. 

Suppose that we have an initial value problem for the Pfaffian system (5-19) (x) and that we consider 
two evolutions, one in which H 1 are constants in the parameters and one in which they are preassigned arbitrary 
analytic functions with values in the prescribed centralizers. Let us denote by ^ and 1 5 the two kernel solutions 
of the Pfaffian system. Define then the function G(z,t) (here t denotes collectively all the isomonodromic 
deformation parameters) as follows 

G(x,t) = *^(x)(^)- 1 (x) . (5-23) 

Since both Stokes' phenomena are the same, this implies that G(z) is a single-valued analytic invertible function 
on the punctured domain CP 1 \ V. Since the essential singularities of * and VP have the same asymptotic 
expansion in the same sectors of the irregular singularities, this implies that G{x) has no singularity there. The 
only possible singularities are poles at the Fuchsian singularities. The same considerations used in Section (4) 
show that in general G(x) has poles of order equal to the resonance index at the given Fuchsian singularity (in 
case of nonresonant singularity G(x) must be analytic there). This shows that the RH problem for G discussed in 
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Sect. 4 admits a solution for arbitrary choice of the group elements H 1 in the centralizer of the local monodromy 
matrices. In other words the dependence of the // 7 's on the parameters t is "pure gauge" and can be completely 
gauged away by means of a rational gauge equivalence G(x,t) which does not alter the singularity structure of 
A(x) . 

In order to fix part of this arbitrariness we can assume that the coefficient A^^^ is in the gauge-fixed 
form (3-92) within each block: this fixes an overall gauge Y°° uniquely (up to the centralizer of ^.oo.r^) 
for our connection. At this point the procedure of Thm. 3.1 for the singularity at oo produces a unique 
formal gauge Yoo in which the leading coefficient is necessarily in the centralizer of Aoo, roa ; the entries 
of this leading coefficient will be constants of the motions exactly as in the previous guide-example, and 
hence could be set to zero in the sense that this would be a consistent reduction of the problem (but of 
no practical advantage). 

Corollary 5.1 (Normalized isomonodromic deformations) Let A(x) be as in Thm. 5.3. Then the 
one form obtained from eq. (5-18) by setting dY°° = = AH 1 defines an integrable Pfaffian system. 

6 Isomonodromic Tau function 

The authors of [11] gave a definition of tau function in terms of a suitable closed differential on the space 
of times. The definition was expressed in terms of (formal) residues of the formal local gauge (which we 
denote Y 1 ) around each singularity of the connection. It is shown in [6] that that definition is equivalent 
to one which manifests the spectral nature of the tau function. In fact the definition in [6] is much more 
convenient and easily generalizablc to the present context 

We first need to fix the gauge arbitrariness of our deformation system by requiring that the deforma- 
tions occur according to the integrable Pfaffian system specified in Corollary 5.1. 

Let us set up the notation: we consider the spectral curve (characteristic polynomial) of the con- 
nection, namely the set of C 2 

n 

E(x,y) := det(y-A(x)) =£y'P n _ J [x;A] = . (6-1) 

3=0 

We think of this curve S in P 1 x P 1 and we regard the two functions x, y : £ — > CP 1 as meromorphic 
functions on the curve, or projections. The x projection x : £ — ► CP 1 has ramification points at all the 

7 S X>i (irregular singularities) and at oo, together with other ramification point whose position is not 
known a priori and in general depends on the " times" . 

If we perform an appropriate desingularization of the curve above each pole x — 7 of the divisor T>\ 
(including 00) we obtain a branching structure in which the n sheets of the projection x : £ — > CP 1 arc 
glued precisely according to the dimensions of the Jordan cells of the leading coefficient of the connection 
and cyclically permuted by a small loop of x around 7. For example, suppose the Jordan form of A(~) )T 
is (we set r 7 = r for brevity since the discussion is local anyway) 

[A {j)>r ] = (A 1 5 ) 2 A 2 3 A 3 (A 4 ) 2 (6-2) 

This means that there are 2 Jordan blocks of size n\ = 5 = n 2 and eigenvalue Ai, one Jordan block 
of size 3 and eigenvalue A 2 etc. (this would correspond ton = 5x2 + 3 + l + lx2 = 16). Under 
our assumptions for the Lidskii matrix of the subleading term, these eigenvalues split into distinct cyclic 
k— tuplets where k is the size of the block and the number of k— tuplets equals the number of blocks of 
that size. In the example we will have two groups of five-plets yij, D2,j, j = 1 . . .n\ = 112 that have the 
common asymptotics 

2/aj ~A 1 (x- 7 )- 1 -'-(l + (l)). (6-3) 
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Both groups admit a Puiseux series expansion in powers of £ := (x — 7) ' ll and each group forms a cyclic 
ni-tuplet 

y a , j (* 2Wna = yaj + i(ti). (6-5) 

The two groups are distinct in the sense that the first coefficients of the expansion in the brackets differ 
from each other c 0) i ^ q>,i, a 7^ b. The first r ■ n a — n a r c coefficients are actually the isomonodromic 
times of our problem, or more precisely (carrying on with this example) 



This means that we can "extract" the isomonodromic times (as well as other parameters of formal 
monodromy) with the residue 



fn a — l 



t a .j = res (u a ) Ja ya.a+i n a £"" ) dx = res y(x - 7) ^ dx 

\<r=0 / ^ 



J = l,...,n a r, u a :=e 2 ™' n " 

n a -l f 71 "' 1 \ 

t a .o = to ^ — = res V" j/a.a+i dx = res ydx (6-7) 

Note that the fractional power is in fact a well-defined local coordinate on the (dcsingularization of the) 
spectral curve. The notation is that is the point on the spectral curve that projects to 7 = x(£") and 
a distinguishes the different (in this case a = 1,2) points projecting to the same 7 (and corresponds to 
different Jordan cells of the same dimension, in this example n a = 5, a = 1,2). 
We now define the isomonodromic tau function by the following equations 

Q -. / T i g - 1 \ 

lnr : = H a , K = - res V y a ^ +1 {uj a ) J %x - 7 ) K /" dx (6-8) 



dt a .K ' K x =1 \ £-~> 

d 1 

— lnr: = - res TrA 2 (x)da; (6-9) 
177 2 x=~y 

Quite clearly the above definition should be repeated for all times and all singularities of our connection, 
where for the Fuchsian singularities the only pertinent equation is (6-9): then we should prove that the 
definition is well-posed, namely that the differential defined by these equation is closed. 

In fact -as it turns out- the proof of this fact is no different from [6] and the previous simple case of 
polynomial connection with single Jordan leading block, needing just some modifications in the notation. 

It is based mainly on the following 

Lemma 6.1 Let A(x) be a rational matrix and let M denote the classical adjoint of a matrix M. Let 
y{x) be the (multivalued) eigenvalue of A(x). Define the following expression 



Tr((y-A)(x) (y> - A) (a/)) 



dxdx' 



^ ■■= / ^ x 7 7 x 72 .a? > ( 6 " 10 ) 



Tr 



(y-A)(xj) Tr((y>-A)(x>)) ( x ~ x ') 2 ' 



where y(x) and y'{x) are two determination of the eigenvalues above the point x. Then Q, is a well defined 
bidiffcrential on the spectral curve, with a double pole without residue on the diagonal away from the 
singularities of A and from the branch-points of y(x). 
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Note that the lemma states that -in spite of the (x — x') 2 denominator, there are no singularities unless 
y(x') — > y'(x) as x' — ► x (i.e. unless y and y' belong to the same sheet of the x-projection). 

Proof. This is essentially the same argument used in the proof of Prop. (3.8). The only -a priori- 
singularities of are at the poles of A and possibly at the branch-points of y -which are of no interest to 
us- and on the points above x = x' . We must prove that if y(x) and y(x') are on different sheets then 
there is no pole and if they belong to the same sheet then there is a double pole without residue. Let x 
and x' be in a common neighborhood which does not include any branch-points; then we can distinguish 
the r sheets of the spectral curve yi, . . . , y r . The expressions 



n Q (x) 



A - y a {x) 
Tr(A^y~ a )(x) 



(6-11) 



are then the rank-one spectral projectors on the (one-dimensional) eigenspace and they are well defined 
in said neighborhood. As we already remarked we have 



Tr(n a (.T)n b (.T')) = s ab + 0((x - x') 2 ) 
which proves the assertion. Q.E.D. 



(6-12) 



Remark 6.1 In fact one may prove a stronger assertion (which is not important for our purposes) that Q is a 
well defined differential on the spectral curve everywhere except the diagonal (x, y) = (x' , y'). 

The bidifferential f2 appears naturally when computing the closure of the differential of the tau function 
and the absence of residue on the diagonal is the key feature that allows proving such closure. 

We now remark that -in the same way as in Section 3.1.1- the formal gauge Y 1 {x){x — j) G ~t W 1 near 
any irregular singularity coincides with the eigenvector matrix up to high order. More precisely, near 
x = 7 (r 7 > 0) we have 



A(x) = {Y^)~ 1 D''(x)Y^ - (yT)- 1 —^ 

da: 



(6-13) 



where D 1 {x) is the block diagonal matrix (we momentarily suppress the explicit reference to 7 from the 
labels for notational convenience) 



/ rn 



D(x) 



J ' 1 ni{x~ 1 ) 



\ 



E 

J=0 



tj. 



Hi 



(G^_ 

x — 7 



,D s (x)) 

1 Tlj 



1 



V 

= diag(L»i(a;), . . . 
G } := diag (0, 

Define now the block-diagonal matrices 

G:=diag(G 1 ,...,G s ) 

W := diag(Wi, . . • , W 2 ) , Wj := [wjV-W-* 



\ 



G, 



x — 7 



i.k 



e j 



(6-14) 

(6-15) 
(6-16) 



(6-17) 
(6-18) 
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Then 

A = Yx- G Wdiag{D 1 {x),... 1 D s (x))W- 1 x G Y- 1 +Y'Y- 1 (6-19) 

^ = ^£ f ^-^ n ^ - ; ( .- 7) (6 - 2o) 

fy: = diag(l,w J - J w i a ,... 1 fa;" , ~ 1 ) (6-21) 

As in Section 3.1.1 this implies that Z := Y(x — j)~ G W coincides with an eigenvector matrix P up to 
order (x — 7) 1 ~ r+e where e = min i^r j an( l m t urn this implies that the deformation matrices are 



JTj. 3 = (ydiag(0 J ...,W/,...,0)y- 1 ) p =((x-7) Zdiag(0 W.il/U, : 0)3 



i>i> 



Pdiag(0, . . . , WjQ/Wj- 1 , . . . 1 Q)P~ 1 J (6-22) 



j = l,...,s, J=l,...,rnj (6-23) 
C = (fflr 1 ) =A 7 (x) (6-24) 



For the deformation of parameters at oo the principal part is to be replaced by the polynomial part and 
then the identity (6-22) is valid only up to the x-independent constant. Let us denote by yj,a{x) the 
eigenvalues which are asymptotic to 

y ^ = n-(x-^ +Wj , \r+i-^ +■■■ > " = 0,. .■,»,•-! (6-25) 

?W 7 J nj ( x - 7) T "3 

Note that these rij eigenvalues are cyclically permuted by a small loop around the singularity 7 on the 
base curve. Moreover suppose now that there are two such cyclic multiplcts with the same order rij and 
same leading term t rn . (which means that the leading term at the singularity of A(x) has two Jordan 
blocks of the same size and the same eigenvalue); then our genericity assumption (2.1) precisely implies 
that the next-to-leading coefficients displayed in (6-25) will distinguish the two cyclic groups. We denote 
accordingly the corresponding spectral projectors 

n jg (g):= {A Zll^ Kx) . (6-26) 
Tr(A-y M )(x) 

Then the previous deformation matrices can be written as (restoring the label 7) 



T h= 7 £(* " ' - I , ^-:=e 2 "/"? (6-27) 




VP 



C 7 = Ay (a:) . (6-28) 

Note that the expressions here above which may in principle contain fractional powers of the local 
parameter, in fact do not because of the cyclicity properties. 

Equation (6-25) (and an analogue for x = 00) could be rephrased in more geometrical terms by saying 
that the desingularized spectral curve £ above each point x = 7 has s(=number of Jordan cells) distinct 
points £y which are all branchpoins of order nj for the map x : X — ► CP 1 , for which a local parameter is 

q-y.j := (x — r )) 1 ^ rLj (or qooj := x^ 1 ^^ ). With these notations in place we can finally prove 
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Proposition 6.1 [Tau function for minimally resonant irregular/Fuchsian isomonodromic deformations] 
The following differential is closed and hence is the differential of a locally defined function 



dinr :=EEr^^^>+EEE^7?rrw + 5 £ s ' ^ 9 ) d ? ( fi - 29 ) 

j=0 J=l 7 eCj=0 J=l V x 7J J ^7677 

Here the residues of the differential ydx are taken at points of the spectral curve (i.e. choosing the 
appropriate eigenvalue y 1M ) and are taken with respect to the local parameters by going around the 
point x — 7 of the base-curve n], times. 

Moreover the isomonodromic times and the parameters of formal monodromy for the irregular singu- 
larities are obtained by the residues 

tffj = resx~ J/n Tydx , j = 1, . . . Soo , J = 1, . . . , n°° (6-30) 

iJ J 

t] . = rcs(.T - 7 ) J ' n lydx , j = 1, . . . , « 7 , J = 1, . . . , n], oo + 7 e V (6-31) 
n 1 - 1 

t< = res ydx +^—, (6-32) 

^7 

Proof. The residue-formulae for the isomonodromic times follow from the expansion of the eigenvalues 
in the respective local parameter. The parameter of formal monodromy to. 7 have a correction which was 
explained in (3-113) and follows simply from 

n 1 - 1 

res TrDj(a;)da; = res ydx =t2 i - - J - — , (6-33) 

7 J Ct 2 

where -Dj >7 was introduced in (6-15). The proof uses the same arguments adopted earlier but it is only 
more involved due to the presence of multiple singularities and times. First of all we note that the residues 
can be pushed down on the base-curve 



1 , 1 n 1 \ 

- res(x - 7) J/n lydx = - res YV:r - 7) J/n huj 1 . k ) Ja y 1 , k , <7 dx 

J cl J x --y \ e=0 ) 



(6-34) 



because the sum in the bracket is a bona fide (local) function of x without fractional powers, due to the 
periodicity of the cyclic nj-plet of eigenvalues. Here the symbol 2/ 7 ,fc 1<T denotes the a-member of the fc-th 
multiplct (corresponding to the fc-th block in the Jordan-cell decomposition) of eigenvalues near x = 7. 
We compute the closure of the differential 

JKd t -y.d t ^ k \a.T = ^res ^ {^ix,h) K<T {x - n) d^y^o J dx = 

" \ <T = / 



X — "f 



J res ( £ ^, fe ^(x - ,i) ^Tr K t , fe , CT (x) ( - [A,T^ 
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res 

X — fJ, 



Uv,H K °{x-n) 



(7 = 




res res 

x=[i z— 7 



E,, Jp, , Kc 

(7=0 p=0 



E 



(z - 7) n i (x - fl) 



■ Tr (II At|fciCT (a;)n 7J - iP (z)) 
(z - a;) 2 



dzdx 



This formula is symmetric because we can exchange the order of the residues; indeed if fi ^ 7 the order 
of residues is certainly irrelevant, whereas if fj. = 7 we can exchange the order because the residue of 
the bidifferential vanishes due to Lemma 6.1. In these computations we have assumed that both 7, \x are 
finite poles, but the proof goes through similarly with minor modification only in the notation and local 
parameter if one or both coincide with the pole at 00 and if one or both the deformations are translations 
of the position of the (finite) pole. Q.E.D. 

As an immediate corollary of the proof we have 

Corollary 6.1 (Hessian of the Tau function) The second derivatives of the Tau function are ex- 
pressed in terms of the spectral curve according to the formula? (for the points at 00 the formula needs 
obvious modifications). 



d 2 lnr 



E 

a=0 p=0 



E 



J P,, 



-^Tr(n MAg (x)n 7J - p (z)) 



(z - xf 



dzdx 



A An example: Airy-like equations 

In this appendix we provide an example of isomonodromic deformation of an irregular singularity which 
is of almost trivial nature but illustrates the machinery developed earlier. We consider an arbitrary 
polynomial V(x) of degree n + 1 and the following scalar ODE (note that the simplest nontrivial case is 
Airy's equation) 



V'(8 x )f = xf(x) 
n+1 

J 



3=1 



The matrix first order ODE associated to this equation is 



(1-1) 
(1-2) 



:= 



-v 2 



'J n+1 



V n +1 



1 

-Vn 
Vn+1 



(1-3) 



The matrix A{x) is in fact the companion matrix of the polynomial V'(y) and thus the spectral curve is 
simply the rational curve (genus zero) 

V'(y) = x . (1-4) 

The only singularity is at x = 00 and it is irregular; the leading coefficient is nilpotcnt with canonical 
form 2 (0) d ~ 2 . In a certain sense this is a non-example of our setting because the Lidskii submatrix of the 
sublcading term does not satisfy our genericity assumption. However it is sufficient to perform a single 
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shearing gauge transformation to recast the problem to a non-resonant one and this feature is the only 
essential ingredient to carry out the analysis as in Sect. 3.1. Following the same steps as in Prop. 3.4 
(with the only difference that now we have to perform the inverse shearing), we obtain that there exists 
a formal solution of the form (q = x~ , t n+ i := n{v n+ i)~ 1 / n ) 



(t n+1 /ny G Y{x)c^x tJ ^ = {t n+1 /n)- G q G Z(q)c T ^W- 1 , G := diag(0, 1, . . . ,n - 1) 



Y(x) := J2 



3=0 



n+l 



Q{x) := £ 



n+l 



n+l 



3=0 



3 = 1 



z := 1 



W 



diagonal matrices. 



-(i-l)(j-i)l 



a; := c 



(1-5) 



Note that the asymptotic representation is unique and is an explicit function of the parameters tj because 
the matrix A is completely determined by them and the diagonal matrices Zj are constructed by a 
determined recursive procedure along the lines of Prop. 3.4. Of the two asymptotic presentations the 
one in fractional powers is the usual one in the context of Airy- like equations, whereas the one in terms 
of integer powers is the one we have used in the present paper. The coefficients tj are the expansion of 
the eigenvalue in fractional powers 



, n+l 



3 = 1 



Sj 



n - 1 



= — res x 

x—oo 



-j/n 



ydx 



res (V'(y))-iyV"(y)dy 

y=oo 



n - j v= 

tn 



res (V'iy)) 1 "dy , j < n 



(1-6) 



Vn+l 

t n+ i = n{v n+1 )~ 1/r ' 



note that to = because the curve is of genus zero, indeed 

n- 1 



to + 



= — res TrDb are (x)dx = — res A(x)dx = 



(1-7) 



where Dt, are is the bare form of the differential equation (see (1-9)); the residues of the traces of Df, are 
and A coincide because they are connected by a formally analytically invertible gauge giving a 0{x~ 2 ) 
extra term which does not contribute to the residue. 

One can write an explicit integral representation of the solutions in terms of Fourier-Laplace integrals 



fk(x):= / e**- v to * jtk (x)= / yi e xv - v W = dif k (x) 



(1-8) 
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where the contours Tk can be chosen in n "homologically" independent ways [4]: in the case n = 2 one 
recognizes the standard integral representation of Airy's functions. The formal asymptotic representation 
in fractional powers as well as the Stokes matrices can be obtained by the steepest descent method, 
whereas the integer powers asymptotics cannot. The bare Pfafhan system is given by 

d x - V ^-W - — | «W = , (d tj - -W (x)) * bare = (1-9) 



3 = 1 



A.l Isomonodromic deformations 



re 



e Q{x) x- (1-10) 



In this case there is no monodromy in the usual sense, only Stokes' matrices, and these are the preserved 
data under the deformation. The parameters of deformations are the tj or -which is the same after a 
change of coordinates- the coefficients vj of the "potential" (this terminology comes from the application 
to random matrices). However our connection A(x) is not gauge-fixed, but it could be done so by 
conjugating it by the constant coefficient of the series Y(x) (which is an explicit function of the tj's). 

The deformation equations are easy to describe because we have an explicit solution. Indeed it is 
immediate from the integral representation (1-8) that 

<>. >i' -<n* (1-H) 

3 

and hence the matrices 

T r .= -±diW-' (1-12) 

trivially satisfy the zero curvature conditions; note that they are polynomials of degree at most 1 except 
for T n+ i which is of degree 2. The deformation equations in terms of the parameters tj can be obtained 
by using the Jacobian of the change of coordinate from Vj to tj. Since these deformations are for a 
connection non gauge- fixed the deformation equations arc in the more general form of Thm. 5.3. 

A. 2 Tau function 

Since the spectral curve is rational, the tau function actually coincides with the tau function of the 
Whitham hierarchy 

n+l 

rllnr = ll .nl , (1-13) 

Hi := -- res x 3/n ydx = - - - res {V'(y)) 1+ ™dy (1-14) 

j x=oo j(n + j) 2/=°o 

It is an exercise left to the reader to check that an integral is 

lnr=-yt j H j n+1 ~ j (1-15) 
2 ^ 3 3 n+ 1 V ' 

3 = 1 

For the case of (translated/dilated/gauged) Airy's equation (n = 2) 

v 3 f" + v 2 f'+v 1 f = xf (1-16) 
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the tau function is 

ln ™ = -^; = 768 W3 4 (1_17) 

In all cases n > 2 it is always a rational expression in the parameters tj or Vj. Note that this equation 
could be transformed to the standard Airy equation by a scalar gauge transformation, a translation and 
a dilation, which would eliminate all the parameters. The first really nontrivial case is the "hyper"-Airy 
equation 

V4.f" + V 3 f" + V2.f + Vlf = Xf , (1-18) 

for which 

1 3^(12^4 + t 2 3 ) 

\nT hAl = . (1-19) 
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